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ABSTRACT 


Spectroscopic and photographic investigations reveal 
a complex azimuthal species structure in the exhaust plume 
of a quasi-steady argon MPD accelerator. Over a wide range 
of operating conditions the injected argon remains col- 
limated in discrete jets which are azimuthally in line with 
the six propellant injector orifices. The regions between 
these argon jets, including the central core of the exhaust 
flow, are occupied by impurities such as carbon, hydrogen 
and oxygen ablated from the Plexiglas back plate of the arc 
chamber. The features of this plume structure are found 
to be dependent on the arc current and mass flow rate. Time- 
resolved spectroscopic velocity measurements, obtained by 
a scanning Pabry-Perot spectrometer system at a current of 
16 kA and a mass flow of 6 g/sec indicate argon ion jet 
velocities of 16,500 m/sec, which considerably exceeds the 
Alfven critical speed of argon. It is found that nearly 
half the observed velocity is attained in an acceleration 
region well downstream of the region of significant electro- 
magnetic interaction. Temperature measurements suggest that 
simple gasdynamic expansion processes in the argon jets can- 
not fully explain the observed acceleration. Recombination 
calculations show that the ionization energy is essentially 
frozen. It is possible that a transfer of momentum from 
the core flow of ablation products to the argon jets may 
contribute to the high argon velocities. 
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CHAPTER 1 
INTRODUCTION 

Electromagnetic plasma accelerators have received con- 
siderable attention in recent years as potential candidates 

(1-4) 

for prime propulsion on advanced deep space missions. 

The advantage of this class of thrusters over other elec- 
tric propulsion devices, such as electrostatic ion engines 

or thermal arc jets, lies in their ability to combine high 

4,2*. 

levels of thrust density (~10 N/m ) with very attractive 

3 

levels of specific impulse (2-5 x 10 sec) . 

Of particular promise is the magnetoplasmadynamic (MPD) 

accelerator in which arc currents typically greater than 

10 A interact with their self-induced magnetic fields to 

exert intense body forces on an injected propellant 
(1-4) 

stream, A typical MPD accelerator is shown m 

schematic form in Fig. 1-1. Under suitable conditions current 
flows between the conical cathode and the coaxial anode in 
a diffuse axisymmetric arc. The interaction of the radial 
component of arc current with the self-induced azimuthal 
magnetic field produces an axial acceleration, which results 
in a "blowing" component of thrust. The axial component of 
the arc current interacts with the magnetic field to produce 
a radial gradient in the plasma pressure which provides a re- 
action force on the cathode surface, giving rise to a "pumping" 
component of thrust. An electrothermal thrust contribution 
also exists, resulting from the unbalanced gasdynamic pressure 
of the resistively heated plasma on the chamber walls. 

Operation of the arc at the highest possible current 
levels is encouraged by the dependence of the electromagnetic 
thrust on the square of the discharge current. ^ However, 
electrode heat transfer problems limit steady-state operation 

3 

to a maximum of about 3 x 10 A. These limitations can be cir- 
cumvented by operating the arc in the so-called quasi-steady 
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mode in which the electrical power is supplied in the 
form of repetitive high current pulses of sufficient dura- 
tion that steady acceleration processes prevail over most 
of each pulse, but which are short enough that the inte- 
grated heat load to the electrodes is tolerable. 

Quasi-steady operation offers additional advantages. 

For example, the MPD operation can be extended to multi- 
megawatt power levels which are completely inaccessible 
to steady-state accelerators and power supplies. Such high 
power operation is of basic interest because it has been 
observed that the overall efficiency of the thruster im- 
proves as the power level is increased . ^ This mode 
of operation also greatly facilitates laboratory study of 
MPD accelerators. The high currents permit the arc cham- 
ber dimensions to be scaled up while maintaining the arc 
intensity in the larger volume, thus considerably easing 
interior access. At the same time, the short duration of 
each current pulse limits energy transfer to modest levels, 
even at the higher discharge powers involved, allowing 
small uncooled diagnostic probes to survive the hostile en- 
vironment in the arc chamber and exhaust plume, and thus 
making possible the difficult task of mapping in detail 
the various plasmadynamic properties over the full interior 
of the discharge. 

Recent experimental investigations of the discharge 
region and exhaust plume of quasi-steady MPD accelerators 
have contributed significantly to the understanding of the 
complex physical processes occurring in this class of thrus- 
ters. Clark, ^ 4 ' 5 ) Oberth^ 6,7 ^ and Turchi^ 8,9 ^ have measured 
the detailed structure of the electric and magnetic fields 
in the acceleration region of the MPD arc at Princeton, and 
established the nature of the current conduction patterns, 
plasma potential profiles and electrode losses. On the same 
device Cory^ 9 ^ has successfully measured exhaust momentum 
flux profiles and chamber pressures, and has determined the 



4 


magnitudes of the electromagnetic and electrothermal com- 
ponents of total thrust. Working with a parallel-plate 
self-field accelerator which retains the essential features 
of the coaxial MPD thruster without flow expansion down- 

( 13 ) 

stream of the acceleration region, Di Capua has studied 
the ratio of energy deposition into kinetic and thermal 
modes and has established that the discharge imparts more 
energy to the flow by resistive heating than by direct body 
force acceleration. His results suggest that electrother- 
mal conversion via the downstream expansion permitted by 
the coaxial geometry of conventional MPD arcs may play a 
significant role in the acceleration process. 

One of the most important indicators of the nature of 
the acceleration mechanisms and of the overall thruster per- 
formance is the exhaust velocity pattern. Consequently 

(2 11-13 30) 

much effort has been applied to its measurement. ' ’ 

Malliaris et al^ 12 " 30 ^ have measured exhaust velocities by 
spectroscopic techniques and have examined the effects of 
arc current and propellant mass flow on the specific im- 
pulse of an MPD accelerator similar to the one investi- 
gated here. The results are instructive but contradictory. 

( 12 ) 

In the earlier work argon exhaust velocities of about 

1.8xl(j m/sec are reported, with the arc operating in the 

so-called "matched" condition. The more recent studies, ^ 3( ^ 

on the other hand, indicate a limiting exhaust velocity 

* (21230) 

equal to the so-called Alfven critical speed ' ' ' - 

3 

8.7 x 10 m/sec for argon - which has been proposed by some 
as a fundamental limit for MPD accelerators. Using elec- 
trostatic time-of-flight probes Jahn and Clark, have 

observed an exhaust velocity of the order of 2.5 xlO 4 m/sec, 
which substantially exceeds both the Alfven critical speed 
of argon .and that computed from the electromagnetic thrust 
relation. ^ A potential shortcoming of the time-of-flight 
method, however, is that it may not record the true stream- 
ing flow velocity because it cannot distinguish between con- 
vective and wave velocity components. If wave propagation 
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components are superimposed on the true convective veloc- 
ity of the flow, the time-of-flight technique will record 
the combined streaming and wave velocities. For the par- 
ticular environment downstream of the accelerator, where 
the magnetic fields are small, the only relevant wave 
system would be the ion-acoustic mode, an electrostatical- 
ly coupled longitudinal oscillation which propagates dis- 
turbances at a hybrid velocity determined by both ion and 

(14) 

electron temperatures. In some regions of the ex- 

haust, where temperatures of the order of 2 eV may exist, 
the wave speed may be commensurate with the true streaming 
velocity, leading to serious overestimates by the time-of- 
flight method. 

These problems, and the lack of general agreement on 
the exhaust velocities produced by MPD accelerators opera- 
ted on argon, clearly point out the urgent need to carry 
out accurate and definitive measurements of velocity by 
means of a technique which neither disturbs the flow nor 
is affected by possible wave propagation phenomena. It 
is the primary purpose of this dissertation to contribute 
to this task. 

To measure the pure streaming velocity components in 
the MPD exhaust it is possible to use spectroscopic Doppler- 
shift techniques. Accordingly, a sophisticated scanning 
Fabry-Perot spectrometer system has been developed, capable 
of recording spectral line shifts and widths on a time scale 
more than an order of magnitude shorter than the length of the 
discharge current pulse. Although this system affords time- 
resolved measurements of both velocity and temperature, with- 
out disturbing the exhaust flow or being subject to wave propa- 
gation effects, it, like any optical technique which relies 
on the self-luminosity of the working medium, is unable to 
provide direct local measurements. This is because observed 
line profiles and shifts result from the integrated effect 
of all emitters along the line of sight and are consequently 
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distorted by spatial distributions of emitter number density, 
temperature and velocity . However, if the exhaust plume 
is azimuthally uniform, local values of velocity and tem- 
perature can in principle he obtained by the application 
of suitable mathematical unfolding procedures , ' 

Previous diagnostic measurements of magnetic fields, 
potential profiles, and impact pressure profiles^ 
have indicated a complete cylindrical symmetry of these 
parameters, and consequently it has been generally assumed 
that the discharge radiance profiles are azimuthally uni- 
form also. 10>15) However, early experiments with the 

(17) 

new Fabry-Perot spectrometer suggested otherwise, ' ' in- 

dicating the existence of a complex but regular azimuthal 
and radial structure in the distribution of the argon pro- 
pellant in the exhaust. The serious implications of this 
apparent lack of azimuthal uniformity immediately motivated 
a thorough photographic and spectroscopic investigation of 
the nature and origins of the species structure of the ex- 
haust plume of the MPD accelerator. The results of this 
study, which has brought to light a number of important 
new features of arc behaviour, are reported in Chapter 3 
of this dissertation. 

The bulk of the present work is devoted to the measure- 
ment of the exhaust velocities and temperatures of singly 

ionized argon, the dominant spectroscopic specie of the 
(15) 

propellant, and to the interpretation of the results. 

The experimental measurements are described in detail in 
Chapter 4, while the various acceleration processes which 
may be responsible for the observed velocities are dis- 
cussed in Chapter 5, Appendices covering ablation phe- 
nomena in the discharge chamber and reviewing the principles 
of the Fabry-Perot spectrometer are included at the end of 
the dissertation. The new information gained in the investi- 
gations described here is significant, providing not only an 
increased understanding of MPD arc processes but also es- 
tablishing important new criteria in the search for optimum 
thruster design. 
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CHAPTER 2 

THE MPD ACCELERATOR FACILITY 
2-1 INTRODUCTION 

The quasi-steady MPD accelerator facility used in the 
present investigation is composed of an arc chamber, a mass 
injection system, a power source, a gas-triggered discharge 
switch, and a large dielectric vacuum tank. A schematic of 
the entire system is shown in Fig. 2-1. In this chapter 
each component is briefly described. Additional details 
are available in the references. 

2-2 THE ARC CHAMBER 

The 12.7 cm diamx 5.1 cm deep cylindrical arc chamber, 
shown in Figs. 2-1 and 2-2, has a conical tungsten cathode 
1.9 cm diamx 2.1 cm long and an aluminum anode with a 10.2 
cm diam exhaust orifice. The inner walls of the chamber 
consist of a Pyrex glass insulator ring and a Plexiglas end- 
plate. Argon propellant enters the chamber at the end wall 
through six nylon injector tubes which are arranged in a 
hexagonal pattern at a radius of 2.54 cm and extend 6 mm 
axially beyond the Plexiglas end wall. The choice of plas- 
tics for the injectors and the end wall is based on their 
ease of fabrication, excellent electrical insulation proper- 
ties, and low initial cost. 

2-3 THE MASS INJECTION SYSTEM 

The propellant mass injection system used in the 

present investigation is one based on valve control, which 

allows operation of the accelerator over a wide range of 

( 18 ) 

mass flow rates. The gas-handling system is shown in 

Fig. 2-1. Fast-acting solenoid valves driven by voltages up 
to 10 times their rated values feed propellant gas pulses to 
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the arc chamber and the discharge switch. Mass flow is 
controlled by the propellant reservoir pressure and by 
calibrated orifices in the valves. For a given orifice 
the propellant feed rate can be varied over nearly an 
order of magnitude. The gas pulse to the switch is de- 
layed through a tube whose length is adjusted to synchro- 
nize switch breakdown with the achievement of steady cold 
gas flow in the arc chamber. 

2-4 THE POWER SOURCE 

The 160 kilojoule power source, which provides a 
nearly square nonreversing current pulse to the arc, con- 
sists of an It— C ladder network of 3300 pfd total capacity 
deployed in 120 stations . ^ 13 ^ (Fig. 2-3) These are as- 
sembled in four lines of 30 identical sections, which can 
be cross-connected in various series-parallel combinations. 
For example, the full series connection can drive currents 
up to 26 kA through a matched load for 1 msec at a maxi- 
mum charge voltage of 8 kV while the full parallel con- 
nection of the four lines can provide 100 kA for 0.25 msec. 

Each capacitor is nested in a plywood pigeonhole in the 
structural frame, which also serves to protect adjacent 
units in case of a capacitor failure. One of the capacitor 
terminals is connected directly to ground, and the other 
to the inductor via a 10 cm length of #16 AWG tinned copper 
wire enclosed in a small bore nylon tube which serves as a 
fuse for that section. Each inductor is a three-turn por- 
tion of a continuous solenoid, 10 cm I.D. x 6.4 cm pitch, 
wound from a 3.8 cm x 0.32 cm copper strip. 

The network is matched to the arc load by an appro- 
priate ballast resistor connected in series with the ac- 
celerator. To accommodate the large energy dissipation 
this resistor is actually an electrolytic cell, consisting 
of copper plates immersed in a cupric sulfate solution. 
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Each of the four network segments can be equipped with its 
own electrolytic resistor, all four of which are needed to 
match the full parallel line configuration „ A typical cur- 
rent waveform delivered by the power source in full series 
configuration, charged to 4 kV, is shown in Fig. 2-6a. 

2-5 THE GAS-TRIGGERED SWITCH 

The power source is connected to the accelerator 

(19) 

through a gas- triggered discharge switch, shown sche- 

matically in Fig c 2-4. It consists essentially of a pair 
of 13.8 cm diam stainless steel disk electrodes separated 
by a Plexiglas insulator ring 12.5 cm I.D. x5.1 cm long. 

In operation the switch is evacuated and the electrode gap 
sustains the full bank voltage until a gas pulse is intro- 
duced through the central manifold. When the pressure 
reaches the Paschen limit at the applied voltage, the switch 
breaks down, closing the circuit to the arc chamber. The 
pressure developed in the switch during the current pulse 
is relieved by an outlet to an evacuated dump tank. 

2-6 THE VACUUM FACILITY 

The accelerator proper is mounted at one end of a 
92 cm I.D. x 183 cm long Plexiglas vacuum tank, shown in 
Fig. 2-5. The dielectric nature of this vessel pre- 
cludes perturbation of the exhaust plume by currents which 
might otherwise flow from it to metallic tank walls. The 
tank is equipped with nine 15.2 cm diam access windows which 
can be used for optical diagnostic experiments, and with an 

internal adjustable platform which can support a variety of 

-5 

diagnostic probes. A 10 torr background pressure environ- 
ment is maintained by a 15.2 cm oil diffusion pump backed 
by two mechanical pumps . 

2-7 TYPICAL ARC CHARACTERISTICS 

Figure 2-6a displays typical voltage and current charac- 
teristics of the arc for a 16 kA x 1 msec current pulse and 
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an argon mass flow of 6 g/sec. After an initial transient 
of about 150 nsec the arc voltage reaches a steady-state 
value of 170V, which it maintains until the current ceases. 
Figure 2-6b displays the quasi-steady current distribution 
at these operating conditions, in the form of enclosed cur- 
rent contours plotted on a cross section of the arc chamber. 
These contours, which represent the fraction of the total 
arc current lying downstream of a given location, were de- 
rived from time integrated magnetic probe signals through- 

( 45 ) 

out the arc chamber. ' ' The current distribution, as well 

( G 7 1 

as floating potential profiles' ' ’ and impact pressure 
profiles show no evidence of current spoking, and ap- 

pear to be azimuthally uniform. 
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CHAPTER 3 

SPECIES STRUCTURE OF THE EXHAUST 
3-1 INTRODUCTION 

One of the most widely used diagnostic tools in the 
study of MPD thrusters has been the visual or photographic 
observation of the exhaust plume. ' ' ' ' ' This is not 

surprising since virtually all types of MPD arcs display 
striking patterns of intense visible luminosity in their 
arc chambers and exhaust plumes. Depending on the propel- 
lant material and flow rate, the prevailing magnetic field 
pattern, and the arc power level, a variety of regions of 
characteristic colors and dimensions are usually found to 
appear with rather sharp definition. Studies of these 
luminous features have played an important role in the de- 
velopment of theoretical models of certain types of MPD 

( 2 , 8 ) 

arcs . 

The quasi-steady MPD discharges under investigation 
here have been extensively photographed with a Kerr-cell 
camera and the observed luminous patterns reported in con- 
siderable details, ^ 8 ' Figure 3-1 displays a Kerr-cell 
photograph^ of the chamber interior at 150 nsec into the 
current pulse, by which time the discharge has achieved full 
quasi-steady operation. A portion of the arc chamber side 
wall is removed, as shown in Fig» 3-la, allowing the entire 
discharge region to be viewed through the glass insulator. 

The shutter remains open for 5 nsec. Color slides corres- 
ponding to the black and white photograph shown here reveal 
several interesting luminous features: 1) a bright green 
layer apparently covering the gas injector ports and most 
of the surface of the Plexiglas insulator, extending about 
one cathode length into the arc chamber; 2) an intense blue- 
white conical region between the cathode and the anode, with 
a sharp upstream leading edge and two horizontal dark reddish 
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bands; and 3) a predominantly bluish plume downstream of 
the anode, exhibiting an interesting radial structure of 
three blue "jets" separated by two dark reddish bands 
which are extensions of the ones in the conical region. 

(4 5) 

Based on the azimuthal uniformity of current, ' po- 
tential^ 6 ' 7 ^ and impact pressure profiles it has general- 

ly been assumed that the luminous patterns are also azimuth- 
(4589 15) 

ally uniform. ' ' ' ' However, close inspection of 

color photographs of the discharge makes it apparent that the 
striated structure of the luminous cone and exhaust plume 
cannot be produced by any reasonable cylindrically symmetric 
distribution of gas emissivities . This seeming lack of cy- 
lindrical uniformity has profound implications for the spec- 
troscopic velocity and temperature measurements, and warrants 
further investigation of the plume structure by more de- 
tailed spectroscopic and photographic techniques. 

3-2 SPECTROSCOPIC OBSERVATIONS 

The highly nonequilibrium conditions usually found in 
all MPD arcs make quantitative spectroscopic studies of the 
plasma properties in these devices very difficult. As a re- 
sult spectroscopic investigations have been used chiefly to 

identify the various molecular, atomic, and ionic species 
. . . . (15) 

participating in the discharge. 7 Although the information 
obtained is primarily of a qualitative nature, considerable 
understanding of the structure of the exhaust plume can be 
gained. 

The spectroscopic studies described here were carried 

out using a three-prism Steinheil G„ H. glass spectrograph 

o 

having a linear dispersion of 5.4 to 40 A/mm over the range 
o o . 

from 4300 A to 6600 A, and an effective speed of f/10. To 

obtain satisfactory exposures of the plates it was necessary 

to expose the film for the entire duration of the discharge* 

The author is grateful to Mr. W„ F„ von Jaskowsky who re- 
corded the spectra and made them available for analysis. 
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(15) 

and the recorded spectra are consequently time integrated. 

However, based on two pieces of evidence of the quasi-steady 

nature of the discharge it can be argued that the lack of 

time-resolution is not a serious drawback: 1) the patterns 

and intensity levels of discharge luminosity recorded in 

Kerr-cell photographs of the discharge remain constant for 

(15) 

almost the entire duration of the current pulse, and 

2) the time-resolved light intensity along any line of sight 
through the discharge is found to be essentially constant 
over more than 70% of the current pulse. The latter is il- 
lustrated in Pig. 3-2 which shows the development of in- 

o 

tensity with time of the 4879.9 A All spectral line on a 
line of sight perpendicular to the discharge axis 4.6 cm 
downstream of the chamber exit. (This record was obtained 
photoelec trie ally using the apparatus described in Chapter 
4.) 

The discharge spectra were recorded with a slit width 
of 2.5 xlO -2 mm on Kodak Royal-X Pan sheet film at several 
locations in the arc chamber and exhaust plume (Fig. 3-3) 
at the nominal operating condition of 16 kA, 6 g/sec. The 
discharge was imaged on the spectrograph with a 19 cm focal 
length, 5.0 cm diam achromatic lens such that the vertical 
image slit extends to a radius of 7.5 cm on both sides of 
the centerline. 

o 

The resulting spectra, covering the range from 4300 A 
o 

to 6600 A are shown in Figs. 3-4 and 3-5. For clarity only 
the more prominent lines have been labeled. The following 
radiating species have been identified: All, C 2 , CII, CIII, 
HI, Nil, Oil and WI. All of these constituents are con- 
tained in the propellant or materials of which the dis- 
charge chamber is constructed. Interestingly, no lines of 
neutral argon could be identified, indicating essentially 
complete ionization of the propellant. However, within 
the arc chamber the radiation from neutral argon may be ob- 
scured by the abundance of lines of other species, particular 
ly those of C 2 in the same region of the spectrum. In fact 
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recent spectra obtained in the near infrared have revealed 
lines of AI in the vicinity of the injector orifices . ^ 1 ) 

The only neutral components detected in the visible light 
spectrograms appear to be HI and C 2 , the latter found exclu- 
sively within the arc chamber, as will be shown later. Re- 
cently, near-ultraviolet spectroscopy^ 22 ^ has revealed the 
presence of AIII in the luminous trumpet-shaped region be- 
tween the cathode and anode and in the exhaust, within one 
anode orifice diameter of the chamber exit. In most of 
the exhaust plume, however. All appears to be the dominant 
propellant specie. 

Ideally, when the arc is properly operated, only the 
argon spectrum should be visible. In fact this is not the 
case, and the abundance of impurity lines of carbon, hydrogen 
and nitrogen indicates that a significant amount of material 
from the Plexiglas* end wall and the nylon** injectors is 
being ingested by the discharge. As will be seen, this has 
profound effects on the structure of the arc. (A detailed 
discussion of ablation phenomena is presented in Appendix A.) 

The spectrogram in Fig. 3-4b taken at a point 3 mm in 
front of the cathode tip exhibits a narrow region of intense 
continuum radiation over the entire visible spectrum, a fea- 
ture not observed elsewhere. This is believed to be bremm- 
strahlung radiation generated within a small region of very 

Plexiglas is a polymer of methyl-methacrylate: ^ 23 ^ 

H CHn 

I I 

- C - C - 
I I 

H C - OCH-, 

II 

0 

l- J n 

* * • / o i ) 

Nylon 66 is a polyamide: v ' 

0 OH H 

II IS I i 

- C - (CH 2 ) 4 - C - N - (CH 2 ) 6 - N - 
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high particle density and temperature which exists off the 
cathode tip. 

The spectrograms in Fig. 3-4c and 3-4d, taken 1.25 cm 
upstream and 2.5 cm downstream of the chamber exit plane re- 
spectively, exhibit an interesting radial intensity profile in the 
spectral lines, which is central to the understanding of the 
structure of the exhaust. This is better illustrated in Fig. 

3-6 which displays an enlarged section of the spectrogram of 

o 

Fig. 3~4d in the 4600 A region. The lines of All exhibit 
three relative maxima of intensity, one on the centerline 
and one on either side at a radius of about 1.5 cm. The lines 
of CII , Oil, and Nil, on the other hand, display two off-center 
maxima at radii which correspond to the minima in the argon 
lines. At the centerline the impurity lines exhibit a shallow 
relative minimum. Furthermore, the impurity lines appear to 
have a smaller radial extent than those of argon, indicating 
a concentration of the impurities in the central core of the 
exhaust. This interesting species inhomogeneity, which can 
be identified with the striated structure of the exhaust lumi- 
nosity in Fig. 3-1 is detectable as far as three anode orifice 
diameters (~30 cm) downstream of the chamber exit. (Fig. 3-5) 

3-3 FILTER PHOTOGRAPHY 

Although the spectrograph is capable of yielding in- 
formation about all the radiating species simultaneously, it 
is severely limited in spatial resolution, allowing examina- 
tion of only one dimension at a time. A further step in ex- 
posing the species structure is to photograph the entire dis- 
charge through narrow spectral filters selected to isolate 
certain lines or groups of lines of the various radiating 
species. By using a number of different lines of sight it 
is possible to obtain a fairly accurate assessment of the 
three dimensional distribution of each radiating constituent. 

Ideally, the filter photography should be done with a 
Kerr-cell shutter to obtain time-resolved information, how- 
ever this is not possible because narrow band filters 
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transmit only a small fraction of the entire visible radiant 
energy (typically 1%) , precluding good exposures on avail- 
able low grain films at microsecond exposure times. 

However, because of the quasi-steady nature of the discharge 
over most of the current pulse, the same arguments which are 
used to justify time-integrated spectroscopy are applicable 
to time-integrated photography. The photographs are recorded 
by simply leaving the camera shutter open for the duration 
of the current pulse. 

The spectral transmission curves of the five optical 

filters used in this investigation are shown in Fig. 3-7. 

All are narrow-band interference filters ranging in half- 
o o 

width from 10 A to 90 A. The prominent spectral lines of 

the discharge which are transmitted by each filter are also 

o o 

shown. The 4880 A and 4486 A filters are used to observe 

o o 

All emission, the 5910 A and 5510 A filters CII and C 0 

o A 

emission, and the 6580 A filter HI and CII emission. 

With each of the filters the discharge was photographed 
at five different lines of sight, three in the horizontal 
plane containing the discharge axis, giving side, perspective 
and end views of the exhaust plume and discharge chamber, 
and two in a plane inclined 30° above the horizontal, also 
containing the axis, giving side and perspective views. 

Figure 3-8 displays photographs of the discharge chamber as 
recorded along these lines of sight. Two cameras were used 
in the study. Black and white photographs on Type 55 P/N 
Polaroid positive-negative film (ASA 50) were obtained by 
means of a 4x5 Graf lex Speed Graphic camera with an f/4„7 
lens. Color slides on High Speed Ektachrome (Daylight, 
ASA160) were obtained with a Nikon F 35mm single lens reflex 
camera equipped with either a 55mm Micro Nikkor Auto f/3.5 
lens or a 200 mm Nikkor Auto f/4 lens . The photographs shown 
in this dissertation are enlargements made from the Polaroid 
negatives. 
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3-4 ARGON DISTRIBUTION 

Photographs of the discharge taken through the 4880 A 
All filter at the nominal operating condition of 16 kA and 
6 g/sec argon flow are shown in Fig. 3-9. It is evident that 
the distribution of All exhibits a complex azimuthal struc- 
ture related to the location of the mass injectors (which 
can be seen faintly in the end views) . Downstream of the 
cathode the ionized argon appears collimated in six discrete 
jets which maintain their identity well downstream of the 
arc chamber exit. (The abrupt transverse cutoff in the 
luminosity about an anode orifice diameter downstream of 

the exit in these photographs is a result of the narrow 
o o 

bandpass (10 A) of the 4880 A filter. The bandpass of an 

interference filter shifts towards shorter wavelengths as 

the angle of the incident light deviates from the filter 
(24) 

normal. If the bandpass is narrow the filter will have 

a narrow acceptance angle and thus a limited field of view.) 

The side and perspective views in the 30° plane (Figs. 
3-9d and 3-9e) show a complete absence of All radiation along 
the axis of the plume suggesting the lack of argon in that 
The central luminous spot visible in the end view. 
Fig. 3— 9c, corresponds to the previously mentioned region of 
intense continuum emission immediately in front of the cath- 
ode tip. The discrete jet structure observed in the photo- 
graphs of Fig. 3-9 explains the curious radial intensity 
structure of the spectral lines of All in Figs. 3-4c - 3-6, 
which were recorded on a line of sight perpendicular to the 
discharge axis in the horizontal plane. The three maxima 
of the spectral lines correspond to the three luminous jets 
visible in Figs. 3-9a. 

Great care must be exercised in interpreting the lumi- 
nous jets to be regions of argon concentration. The emis- 
sivity of a given radiating specie for a particular spectral 
transition is not only a function of its number density, but 

also a strong function of the prevailing electron tempera— 

(25 26 ) 

ture. ' Thus the observed patterns of luminosity could 
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conceivably be a convolution of the electron temperature 
field and the species density distribution. Fortunately, 
measurements of electron temperature ^ indicate that it 
has a uniform value of about 1.5 eV throughout the 
entire arc chamber and the near exhaust. Furthermore, in 
the exhaust plume there appears to be only a gentle nega- 
tive axial gradient in electron temperature, as will be dis- 
cussed in Chap. 4. Consequently, it can be assumed that 
the observed luminous patterns correspond to maps of species 
density. 

o 

Photographs through the 4486 A filter (Fig. 3-10) reveal 

a luminous structure virtually identical to that observed at 
o 

4880 A, a result which is consistent with the similarity of 

the radial intensity profiles of all the spectral lines of 
o o 

All. The 4486 A filter is broader than the 4880 A filter 
o o 

(90 A vs. 10 A) and thus transmits more continuum light and 
a few impurity lines (Oil) causing a slight loss of definition 
of the jets. 

Recent spectroscopic investigations in the near-ultra- 
( 22 ) 

violet reveal a distribution of AIII similar to that of 

All, providing compelling additional evidence that the All 
jets do actually define regions of argon concentration. With 
the apparent complete absence of AI from the exhaust plume 
it must be concluded that there is in fact very little or no 
argon between the jets and in the central core of the exhaust. 

Flow Visualization 

It is important for the velocity measurements discussed 
in Chap. 4 to determine whether the observed argon jet struc- 
ture is bounded by flow streamlines, i . e : whether the flow 
direction is along the luminous jets. A simple test is to 
place a spherical body in one of the argon jets and note the 
direction of the resulting bow shock* To a first approxima- 
tion this will reveal the local flow direction. Figure 3-11 

It is known from a previous series of wedge flow studies v 
that the flow in the exhaust plume is everywhere supersonic. 
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displays the results of such an experiment, photographed at 
4880 A. An 8 mm O.D. spherical bulb blown on the end of a 
5 mm O.D. pyrex glass tube is immersed in the upper near 
argon jet in Fig. 3-lla producing a distinct bow shock. The 
axis of symmetry of the shock appears to be parallel to the 
luminosity, indicating that the local flow direction is along 
the luminous jet. Figure 3-llb shows the probe in the upper 
far jet, with similar results. If the probe is placed on the 
exhaust centerline. Fig. 3-llc, no shock is visible through 
the argon filter, attesting to the lack of argon in the core 
of the exhaust. When the probe is moved horizontally into 
the middle near jet. Fig. 3-lld, the bow shock reappears. 

Since the axis of symmetry of the shock wave is in this case 
parallel to the discharge axis, it implies that there is no 
rotational component of velocity in the exhaust plume. 

It is interesting to note, at this point, that at the cham- 
ber exit the radius which encloses the six luminous argon jets 
also encloses virtually all of the flow momentum exiting from 
the chamber, It is not clear whether all the injected 
argon is also contained within this radius. It has been sug- 
gested by Cory^ 10 ^ that as much as 40% of the injected pro- 
pellant may be leaking out of the chamber at low velocity 
through the annular region between the luminous plume and the 
anode orifice. That no significant argon radiation from any 
level of ionization can be detected spectroscopically in this 
region would tend to cast some doubt on this contention. 

3-5 DISTRIBUTION OF IMPURITIES 

In Fig. 3-6 it can be seen clearly that the lines of CII, 
Oil, and Nil have a smaller radial extent than those of All 
and exhibit maxima where the lines of All exhibit minima. It 
can be concluded from this that the regions devoid of argon 
must be occupied by these impurities which are produced by the 
ablation of the Plexiglas end wall of the arc chamber. Again, 
by photography through suitable filters, the distribution of 
the ablated material in the exhaust can be determined. 



37 


Distribution of Carbon 

Figure 3-12 displays photographs taken through the 
o o 

5910 A, 75 A bandwidth interference filter at the 16 kA, 

6 /sec operating condition. Prominent in this spectral band 

o o 

are the 5889.97 A and 5891.65 A lines of CII and the tail of 

o 

the C 2 0-2 Swan band system whose head is at 6191.2 A. It 
can be seen, by comparison with the photographs in Fig. 3-9 
and 3-10 that these constituents are found precisely in those 
regions of the chamber and exhaust where argon is absent and 
vice-versa. The luminous patterns inside the chamber are 
generated by both C 2 and CII, while those outside in the ex- 
haust plume are generated by CII alone, as will be shown 
later. The side and perspective views exhibit clearly the 
expected concentration of the CII near the exhaust centerline. 
The slight decrease in luminosity at the centerline visible 
in Figs. 3-12a and 3-12b was noted previously on the spectral 
lines of CII and other impurities. When viewed along the 30° 
inclined plane, which cuts between injectors, the CII plume 
shows a distinct increase in luminosity at the centerline. 

These interesting features are a result of looking at the 
star-shaped cross section of the CII plume on edge from the 
two different sets of viewing angles. If the plume is viewed 
along the horizontal plane the line of sight across the axis 
sees a smaller thickness of CII than on either side. On the 
other hand, at an angle 30° to the horizontal the line of sight 
across the axis passes through two diametrically opposed arms 
of the star shaped cross section of the CII flow and the center- 
line appears brighter. 

Flow visualization (Fig. 3-13) by means of the spherical 

probe technique used in the study of the argon jets helps 

confirm some of the conclusions drawn above. Figure 3-13b„ 

o 

photographed at 5910 A, shows the 8 mm spherical probe im- 
mersed in the carbon plume on the centerline of the exhaust, 

4.6 cm from the anode plane. A very strong bow shock is 
produced by the flow over the probe, confirming the 

presence of carbon in the core of the exhaust. (The shock 
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angle corresponds roughly to a flow Mach number of 2.2 . 

If it is assumed that the exhaust flow can be approximated 

by a source flow model the flow Mach number is estimated 

to be about 6 where the shock reaches the visible boundary 

of the carbon jet.) When the probe is located in one 

of the argon jets, as shown in Fig. 3-13a, no shock is 

o 

visible at 5910 A, indicating the absence of carbon there. 

The star-shaped luminous pattern visible inside the 
chamber in the end view of the discharge clearly indicates 
where ablation of the Plexiglas back plate occurs. The 
ablation pattern etched on its surface can be easily dis- 
cerned in the end view of the chamber in Fig. 3-8. In the 
perspective views in Fig. 3-12 it can be seen that a luminous 
star-shaped layer extends about a cathode length into the 
chamber from the back plate. This immediately identifies 
it with the luminous layer along the end wall visible in the 
side-on Kerr-cell photograph of Fig. 3-1. The products of 
ablation generated at the insulator surface over this star- 
shaped region retain this geometry as they are accelerated 
downstream and are thus responsible for the six- jet struc- 
ture of the argon in the plume. 

In Fig. 3-14 are shown three enlarged sections of the 

spectrograms of Figs. 3-4a, c and d in the range of 5300 A - 

o 

5900 A, which contains the 0-1 Swan band system of . The 
first spectrogram, obtained 1.1 cm in front of the Plexiglas 
back plate of the arc chamber, i . e ; within the luminous star- 
shaped layer at the plate, shows that the lines of C 2 , CII, and 
CIII have similar radial intensity profiles, which indicate 
that these three species have similar spatial distributions 
there. No lines of neutral carbon could be identified, a 
consequence of the relatively small number of lines of Cl in 
the visible spectrum and the presence of many other lines 
which make unambiguous identification very difficult. It seems 
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reasonable however, that neutral atomic carbon must be present 
in any region that contains both the molecular and singly 
ionized specie. 

The rather abrupt termination of the luminosity of the 
layer at the Plexiglas wall about one cathode length from the 
wall is thought to be due to the rapid disappearance of C 2 
through dissociation. Except for a thin layer along the sur- 
face of the circular pyrex wall there appears to be no molecu- 
lar carbon present downstream of the cathode. In part this 

can be seen in Fig. 3-15 which shows two perspective views 
, o 

and one end view of the chamber interior through the 5510 A, 
o 

80 A wide interference filter, which transmits part of the 
0-1 band system. Only a short plume extending from the cath- 
ode to the plane of the anode is visible. This short plume 
does not itself contain any C 2 as can be ascertained from the 
spectrogram in Fig. 3-14b, obtained 1.25 cm upstream of the 
exit plane. All the spectral lines except those of C 2 show 
a marked increase in intensity near the centerline. The in- 
tensity of the C 2 bands remains nearly constant across the 
chamber diameter indicating that the C 2 radiation originates 
either behind or in front of the plume, most likely at the 
cold chamber walls where carbon deposits have been found. The 

short plume in Figs. 3-15a and 3-15b must therefore consist of 

o 

other species which radiate within the bandpass of the 5510 A 

filter. None of these constituents have been positively 

identified but they are believed to be impurities in the higher 

stages of ionization. The spectrogram of Fig. 3-14c shows a 

complete absence of molecular carbon bands 2.5 cm outside the 

o 

chamber. The 5696.0 A line of CIII is also absent, while the 
o o 

5889.97 A and 5891.65 A lines of CII are prominent. Spectro- 
grams farther downstream (Fig. 3-5) exhibit only the singly 
ionized specie of carbon. 

Distribution of Oxygen and Nitrogen 

No photographs were taken to directly reveal the distri- 
bution of Oil or Nil because virtually all of their line radi- 
ation is emitted too close to lines of CII or All. This would 
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o 

require the use of very narrow bandwidth filters ( < 2 A) which 
would severely limit the field of view due to the very small 
angle of acceptance of such filters. However, Oil and Nil 
are expected to have concentration profiles similar to CII be- 
cause their spectral lines exhibit the same radial intensity 
features. (Fig. 3-6) 

Distribution of Hydrogen 

In all the spectrograms of the discharge (Figs. 3-4 and 
o 

3-5) the 4861.33 A H« line of atomic hydrogen is conspicuous. 

, . “ o 

The Hot line at 6565.82 A falls near the limits of the spectral 

sensitivity of the Royal-X Pan emulsion in the red and is thus 

poorly exposed. The Hy line at 4340 A is intrinsically fainter 

than the other two lines ' and apparently so broad that it 

is not easily recorded either. Figure 3-16 shows photographs 

o o 

of the discharge taken through the 6580 A, 90 A half-width 

interference filter which transmits the 6562.82 A H« line 

o o 

and the 6578.03 A and 6582.85 A lines of CII (Fig. 3-8). 

Arc current and propellant mass flow are 16 kA and 6 g/sec 
respectively. It has already been demonstrated that CII is 
found predominantly in the core of the exhaust flow. Thus in 
Fig. 3-16 the high luminosity of the core of the exhaust is 
due in large part to the CII radiation. The fainter but more 
extended part of the plume is due to H * radiation exclusive- 
ly. The hydrogen appears to be present throughout the entire 
exhaust plume, even within the argon jets, and is found at 
radii beyond the visible outer boundaries of the argon and 
carbon flows. Furthermore, inside the arc chamber the distri- 
bution of hydrogen appears to be more uniform than that of 
carbon. This is probably a result of the greater mobility of 
hydrogen due to its very low atomic weight. 

Since ionized hydrogen does not have a line spectrum, its 
presence must be inferred from that of the neutral specie and 
the fact that the other constituents of the exhaust (all of 
which have approximately the same ionization potential as hy- 
drogen) are present in only singly or doubly ionized form. 



FIGURE 3-16 

AP25-P-455 



SIDE VIEW 


PERSPECTIVE VIEW 


END VIEW 

I- 8927 


1-8936 


1-8934 


1-8927 


Reproduced from 
best available copy 


1-8932 


T. QQ9Q 


DISCHARGE AT 6580 A SHOWING DISTRIBUTION OF HI AND Cl. J = 16 kA , m = 6g/sec 


Ui 






46 


The fact that neutral atomic hydrogen is detectable, whereas 
other neutral atomic species are not is most likely a conse- 
quence of the intrinsically greater strength of the atomic 
hydrogen spectral lines. Very small concentrations of 
neutral hydrogen can give rise to the observed intensities. 

The pervasive presence of this impurity is on one hand use- 
ful, since it allows determination of electron densities in 
the chamber and exhaust through measurements of the Stark 
broadening of its spectral lines. A typical contour map of 
electron density obtained in this fashion is shown in Pig. 

3-17 . It is interesting to note that the sharp upstream 

edge of the conical plume between cathode and anode visible in 
all the discharge photographs lies, approximately parallel to 
the local electron density contours, lending evidence to the 
hypothesis of Turchi v ' ' that a trailing shock off the cathode 
tip may be responsible for the sharp upstream plume boundary 

A disadvantage of the hydrogen contamination is that 
even relatively small concentrations of the ionized specie 
can significantly affect every process that involves the dif- 
fusion of positive charges in the plasma. For example, 
measurements of flow velocity by means of biased double electro- 
static probes aligned parallel and perpendicular to the flow are 
sensitive to the concentration and atomic weight of the ions 
impinging on the probe surfaces. If the probes are im- 
mersed in„ a flowing mixture of ions of different mass it is 
necessary to know their relative concentrations to enable re- 
duction of the observed ion saturation currents to flow vel- 
ocities. Another effect of the presence of light ions such 

as hydrogen is a significant increase in the ion acoustic 
( 29 ) 

velocity, which would increase the possible error in time 

of flight velocity measurements ^ 11 ^ if wave propagation effects 
are present. 
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3-6 OPERATION AT OTHER CONDITIONS 

Various investigators have examined the behaviour of 
the quasi-steady MPD accelerator over wide ranges of input 
mass flow, m, and arc current, j # ( 4-7 * 10-13 # 30) R esu it s sug- 
gest that a number of physical processes and terminal proper- 
ties are functions of the parameter J 2 /m, rather than J and 
m separately. Malliaris et al, have identified the so- 

called "matched" operating condition, J 2 /m « 40 kA 2 sec/g, at 
which thruster performance is optimized. (The data presented so 
far in this dissertation were obtained at this condition.) 

Their results indicate that for the "starved" operating con- 
dition, j /m >40, insulator and electrode erosion is a prob- 
lem, while at the "overfed" condition, J 2 /m<40, some of the 
input mass bypasses the acceleration region. It seems reason- 
able, then, to expect the species structure of the exhaust plume 
to be dependent on the operating condition and be character- 
ized by the quantity J /m. To explore this possibility the 
arc was operated over a broad range of current and mass flow, 
and the resulting argon and ablated carbon exhaust flows 
photographed through the appropriate spectral filters. A 

effect on the species distribution has been confirmed. 
Some of the results are presented in Figs. 3-18 - 3-21 in the 
form of perspective and end views, at five different operating 

conditions: J /m s* 10, 20 , 40 , 80 , 160 . The matched condition, 

2 • w 
J /m s® 40 is included for reference. Figures 3-18 and 3-19 

show that if the arc is overfed the argon jet structure be- 
comes less distinct, each jet spreading azimuthally and to- 
wards the plume axis. At the most overfed condition shown 
(Fig. 3-18a and 3-19a) the discrete jet structure is no longer 
detectable, and the argon distribution is much more uniform 
than at the other conditions. When the discharge is starved 
the argon jet structure is clearly retained. At the most un- 
derfed condition shown (Figs. 3-18e and 3-19e) there appears 
to be some loss of definition of the jets, possible due to an 
increase in the continuum emission throughout the exhaust 
plume . 
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o 

The photographs at 5910 A, of ablated carbon in the dis- 
charge (Figs. 3—20 and 3-21)/ show complementary changes. 
Operation at overfed conditions decreases the amount of car- 
bon in the exhaust indicating a reduction in the ablation 
rate. At starved conditions ablation is severe and carbon 
is found in abundance throughout the flow field. At any con- 
dition ablation appears most severe in the annular region 
between the cathode and the injectors. Bright narrow streaks 
visible in some of the perspective views are due to glowing 
particles ejected from the cathode. 

In general, based on these and numerous other photo- 
graphs, it can be said that under all but highly overfed con- 
ditions, i , e : for J /m>10, the exhaust plume of the MPD arc 
exhibits varying degrees of species separation. The tendency 
is for the species structure to be most sharply defined at 
matched or somewhat starved conditions. The occurrence of 
localized ablation of the Plexiglas back plate under even 
overfed conditions indicates that the overall operation of 
the arc cannot be strictly labeled "overfed" , "matched" , or 
"starved". Under a wide range of conditions the arc is ap- 
parently always starved in the regions between the injectors 
and between the cathode and the injectors. In addition, it 

does not appear possible to characterize all features of the 

2 

plume structure solely by J /m over the entire possible range 
of this parameter. For example, the distribution of abla- 
tion products appears to be a weaker function of m at con- 
stant J than might be expected from a J /m scaling re- 

lationship. At highly overfed ( J /m < 10) and highly underfed 

2 • 

(J /m > 160) conditions the input mass flow appears to have no 
detectable effect on the rate of ablation or the impurity dis- 
tribution in the exhaust. A further discussion of ablation 
phenomena is presented in Appendix A. 
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CHAPTER 4 

VELOCITY AND TEMPERATURE MEASUREMENTS 
4-1 INTRODUCTION 

The complex radial and azimuthal species structure of 
the exhaust flow of the MPD arc gives rise to difficulties 
in the interpretation of velocity measurements obtained by 
any of the usual probe techniques. In particular, it raises 
additional questions about the previously published time-of- 
flight velocity measurements^'*'^ where it was assumed that 
only argon existed in the exhaust so that velocities measured 
along the plume axis could be ascribed to this specie. As 
was shown in the previous chapter for all but very highly 
overfed conditions (which are not of propulsion interest) 
there is no argon on the plume axis, and consequently the 
reported time-of-f light centerline velocities correspond 
to the velocities of the ablation products, rather than of 
the propellant. Additional complications in the interpre- 
tation of time-of-flight velocity data arise if ion acous- 
tic wave propagation effects are important. The species 
inhomogeneity would lead to radial profiles of wave veloc- 
ity which are sharply peaked on the centerline because of 
the significantly lower atomic weight of the flow there. 
Unfolding of wave and convective velocities would thus be 
considerably complicated. These problems, typical of 
measurement techniques which cannot distinguish between 
atomic species, give added impetus to the application of 
spectroscopic Doppler-shift techniques to measure propel- 
lant velocities. 

Spectroscopic methods have been used by a number of 
other researchers in their investigation of various types 
of steady self-field and applied field MPD arcs . ^ 16 ’ 31-37 ) 
Most of these devices were operated in low power regimes. 
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with propellants other than argon, and in certain cases ' 
under relatively high environmental pressures. In all but 
two of these investigations^ 6 ' 3 '*'^ no attempt was made to 
reduce the spectroscopic data to local flow properties. The 
results of these studies, although interesting in their own 
right, are only of academic interest in the present context. 
Spectroscopic measurements on quasi-steady high power MPD 
arcs have also been made. ( 12 ' 18 ' 30) Of particular interest 
for comparison with the present work are the results of 
Malliaris et al,^ 12 ' 30 ^ obtained on an accelerator quite 
similar to the one investigated here. 

This chapter is concerned with time-resolved spectro- 
scopic velocity and temperature measurements of singly 
ionized argon, the dominant radiative propellant specie in 
the exhaust of the quasi-steady MPD arc. Unfortunately, the 
lack of cylindrical symmetry of the propellant in the exhaust 

flow precludes the reduction of Doppler-shift measurements to 

(16, 25,38) 

purely local values by Abel inversion techniques. 

However, it will be seen that the rather sharp collimation of 
the argon jets does allow deduction of reliable values of 
average velocity and temperature as functions of axial dis- 
tance within the individual argon jets. 


4-2 EXPERIMENTAL APPARATUS 

Constraints on Spectrometer Resolution 

The resolving power of the spectrometer used in the 
velocity and temperature measurements must be sufficient to 
allow clear separation of Doppler shifted line components 
and also be such that the instrumental broadening is less 
than the expected Doppler line widths . The order of magni- 
tude of the expected line shifts and widths is estimated be- 
low. 

If an emitter of radiation moves with a velocity v 
relative to an observer, then a wavelength X will be observed 
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such that ^ 39 ^ 



X 0 I -p cose 


where Z 3 = v/c, c is the speed of light, 0 is the angle 
between the velocity vector and the line of sight, and A 
is the wavelength which would be observed if the emitter 
were at rest. For non- relativistic velocities, i . e . , 
if « l, the following approximation is valid: 

aX = ht cos 9 

c 

where aX = X 0 — X « X is the Doppler shift. Velocities 
of the order of 10 4 m/sec, which are characteristic of MPD 
devices, lead to Doppler shifts of about 0.1 -0.2 A in the 
visible spectrum. The Doppler broadening, aA d , of spec- 
tral lines due to the random thermal motions of the radiating 
particles is given by^ 40 ^ 

AX D - 7/6 x 10 X 

where A is the wavelength at line center of the spectral 
line, and T and M are respectively the temperature (°K) 
and molecular weight of the radiating particle. (A more de- 
tailed discussion of this relation is given in section 4-4.) 

4 o 

Argon temperatures of the order of 10 K ( ~ leV) , which are 

usually found in plasma accelerators, lead to line widths of 

o 

the order of 0.05 A. 

The high resolving power necessary to detect Doppler 
shifts and line widths of the magnitudes estimated above can 
be achieved with a Fabry-Perot interferometer. ^ 24 ' 41-45 ) 

A review of the theory of this device and of the factors in- 
fluencing its performance is presented in Appendix B for the 
benefit of the reader who may not be familiar with its princi- 
ples. Additional details can be found in the literature . ^ 24 # 41 “ 45 ) 
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The Fabry-Perot interferometer is not only capable of much 
higher resolution than a prism or grating spectrograph* 
but also has a much higher luminosity „ Furthermore* it is 
simple in principle, readily adaptable to time-resolved 
operation and is relatively inexpensive . 

The Fabry-Perot Spectrometer System 

The interferometer and supporting optics used in the 
present study are shown in schematic and photograph form 
in Figs . 4-1 - 4-3 0 The Fabry-Perot (henceforth denoted 
F.P.) is a commercial unit manufactured by Coherent Optics, 
Inc. with interferometer plates by Perkin-Elmer Corp. The 
plates, 38 mm diam x 9.5 mm thick are made of fused silica 
and are polished flat to X /200 ( A « 5000 A) . Their re- 
flecting coatings, which consist of multilayer dielectric 

films, have a reflectivity of 95 ± 1% in the spectral range 

o 

of 4600 - 5400 A. These coatings reduce the effective sur- 
face flatness to about X/lOO . (The reflectivity finesse 

F R and the flatness limited finesse F^ (both defined in Appen- 
dix B, section B-3) are thus approximately equal (F R »F D »50) 

so that the effective instrumental finesse F and trans- 

e 

mission factor °C are optimized with respect to each other.) 
One of the interferometer plates is fastened to one end of 
a hollow tubular piezoelectric transducer. Application of 
a linear ramp voltage to the transducer causes it to trans- 
late the plate perpendicular to its surface so that the 
spacing between the interferometer plates changes linearly 
with time. This effects a linear change in the transmission 
wavelength of the F.P. causing it to act as a tuneable wave- 
length filter which scans a region of the spectrum. (Appendix 
B, section B-2) „ The other stationary plate of the F.P. is 
held in a gimballed mount which can be adjusted by means of 
micrometer screws and piezoelectric transducers to render the 
plates parallel. 

The use of the F.P. as a spectrometer requires supporting 
optical elements and electronic equipment to collimate light 
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from the plasma under investigation and detect the F.P, 
output (Figs. 4-1 - 4-3) . Light from the discharge is 
directed by the mirror M^ through an adjustable diaphragm 
to the objective lens (50 mm diam x 1270 mm f.l.)whose 
focal point lies on the discharge axis. Mirror M-^ can be 
adjusted about both horizontal and vertical axes by means 
of micrometer screws to obtain the desired line-of-sight 
through the exhaust plume. Lenses L 2 (50 mm diam x 190 mm f.X.) 
an< 3 L^ (18 mm diam x -47 mm f.l.) reduce the collimated in- 
coming beam diameter by a factor of four so that only the 
central portion (~1 cm diam) of the F.P. plates is used. 

(This tends to minimize the effects of plate distortion 
caused by accelerations of the moving plate during a scan.) 
Mirror M^ directs the output of the interferometer to the 
lens L^ (50 mm diam x 172 mm f.l.) which focuses the center 
of the scanning fringe pattern on the 0.8 mm diam pinhole , 

PH. The size of this aperture is chosen so that its finesse 
F p matches the plate flatness finesse F Q and the reflectiv- 
ity finesse, F R , as explained in Appendix B fl Lens L g 
(12.5 mm diam x 65 mm f.l.) collimates the light from the 
pinhole for passage through a 4880 A, 10 A bandwidth inter- 
ference filter, IF, which selects the 4880 A line of All for 
investigation. An RCA 1P28 photomultiplier tube detects the 
scanned spectrum. A TRW Model 71-BR pulsed argon ion laser 
(pulse width: 40 psec) provides the reference wavelength at 
4879.9 A against which Doppler shifts are measured. (The 
reference wavelength cannot be obtained from the discharge 
at right angles to the flow velocity, as is usually done, for 
reasons which will become clear later.) The laser output is 
expanded and collimated by lens system L c -L_ (5 mm diam x - 4 mm 
f.l. and 19 mm diam x 67 mm f.l.) and superimposed on the 
collimated discharge light by means of mirror M 2 and a clear 
flat glass beam splitter, BS . The expanded laser beam 
diameter is adjusted by means of an iris diaphragm, ID 3 , to 
match the diameter of the discharge light beam. A set of 
polarizers and neutral density filters is used to control the 
laser intensity reaching the photodetector. 
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All the components of this system are mounted on an opti- 
cal bench made from heavy aluminum channel stock which is 
securely bolted to a 2.5 cm thick x 152 cm long x 61 cm wide 
aluminum plate. The latter is placed atop a 10 cm thick 
plywood slab of the same length and width which is covered 
with a 12.5 mm thick layer of fiberglass vibration isolation 
matting. The wooden slab is cradled in a steel frame which 
in turn rests atop a heavy angle-steel framework. Leveling 
is accomplished by means of four large diameter screws which 
support the frame around the wood base. The entire assembly 
sits atop four 10 cm x 10 cm x 5 cm fiberglass vibration 
isolating pads, which are very effective in preventing floor 
vibrations from reaching the interferometer. The optical 
system atop this optical table is enclosed in a 1.5 mm thick 
aluminum box, lined with a 2.5 cm layer of black foam rubber, 
which keeps air currents, sound vibrations and extraneous 
light from the interferometer and other components. 

Method of Operation 

Proper synchronization of the F.P. scan and laser refer- 
ence pulse with the discharge current pulse is a crucial ele- 
ment in the operation of the spectrometer system described 
above. A schematic of the electronic synchronization and re- 
cording system is shown in Fig. 4-4. The sequence of events 
is sketched in Fig. 4-5. Initiation of the discharge is 
sensed by a Rogowski coil ^ embedded in the anode re- 
turn cylinder. The resulting signal, triggers time base 

A of the Tektronix 555 oscilloscope. This time base, which 
records the discharge current history, completes its sweep 
in 1 msec. The sawtooth voltage (0-150 V in 1 msec) used in 
the horizontal sweep of this time base is applied to the 
piezoelectric translator attached to the moving mirror of the 
F.P., which as a result moves a distance of about 1.7 p to- 
wards the stationary mirror, for a scan of about three free 
spectral ranges in the middle of the visible spectrum. Ideally 
the moving mirror should travel at constant velocity from the 
beginning of the F.P. scan at discharge initiation to the end 
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of the ramp voltage 1 msec later, so that the plate separation, 
and thus the transmitted wavelength is a linear function of 
time. However, the piezoelectric transducer and the moving 
plate have finite dimensions and mass so that the plate motion 
is not linear; rather it has superimposed on it an oscillatory- 
component as shown in Fig. 4-5b. This is tantamount to a loss 
of effective F.P. finesse. Fortunately the oscillations damp 
out quickly, so that during the second half of the scan the 
plate motion is essentially linear in time. Another result of 
the initial acceleration is a bending or bowing of the moving 
plate which seriously degrades the plate flatness and thus the 
overall finesse of the F.P. during the first 100 psec of the 
scan. Some bowing appears to exist throughout the scan, be- 
cause the experimentally measured overall finesse is always 
somewhat less when the instrument is scanned at the present 
rate of 1 scan/msec, than if it is scanned one or two orders 
of magnitude more slowly. Results indicate that during the 
second half of the 1 msec scan the plate flatness is degraded 
to approximately A /90 from the stationary value of A /100, 
which is quite tolerable. 

Superimposed on the ramp driving voltage is a variable 
D.C. bias obtained from a battery operated power supply. The 
purpose of this bias is to provide a fine adjustment of the 
initial F.P. plate separation and thus permit proper centering 
of the transmission peaks on the cathode ray tube of the os- 
cilloscope. The wavelength scan detected by the photomulti- 
plier tube is exhibited on time base B of the oscilloscope, 
swept at twice the rate of time base A ( i. e ; in 0.5 msec) but 
delayed with respect to it by 0.45 msec. In this manner only 
the second half of the F.P. scan is displayed. The D.C. bias 
to the F.P. is adjusted so that the spectral line is swept out 

twice in one sweep of time base B. (Fig. 4-5d) . With the free 

o 

spectral range of 0.78 A used in the present work a time reso- 
o o 

lution of 0.003A/psec at about 5000 A results. As will be 
seen this corresponds to the scan of a typical spectral line 
in about 50 (jisec. 
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The laser used as a wavelength reference can produce only 
a single pulse 40 psec wide during the 1 msec scan of the inter- 
ferometer: thus, it must be triggered at a suitable point in the 

scan such that the plate separation is correct for transmission 

o 

of the laser wavelength (487 9.9 A), as illustrated in Fig. 4-5d. 
Since only the second half of the F.P. scan is recorded, the 
laser pulse must be delayed at least 450 psec with respect to 
discharge initiation. The proper delay is obtained by means 
of a Hewlett-Packard Model 214A pulse delay generator which is 
triggered by the initiation of the sweep of time base A on the 
oscilloscope. 

The overall finesse of the Fabry-Perot spectrometer system, 
when operated as described above, has been measured to be 
F t = 25. Stability of the instrument against external influences 
is very good, except when the room temperature is not constant. 
Before any data are recorded the various pieces of electronic 
equipment are allowed to warm up for at least two hours, or 
until the interferometer displays a drift stability of more 
than 5 minutes. On particularly good days stabilities as long 
as 30 min. have been observed. Clearly, a constant temperature 
enclosure for the interferometer would greatly increase its 
stability. 

4-3 VELOCITY AND INTENSITY MEASUREMENTS 

By proper choice of lines of sight through the exhaust 
plume of the MPD accelerator, a variety of quantitative mea- 
surements may be carried out by means of the system described 
in the previous section. The spatial location of the argon jets 
can be accurately mapped by recording line intensities as a 
function of distance both parallel and perpendicular to the 
discharge axis. Doppler-shift velocity measurements obtained 
at two sets of lines of sight at two different angles to the 
discharge axis can yield information not only on the develop- 
ment of average velocity along each jet, but also on the vari- 
ation of average jet flow angle. Finally, argon ion tempera- 
tures can be obtained from the Doppler broadening of spectral 
lines. All measurements reported here are confined to the 
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"matched" operating condition of J = 16kA, and m=6 g/sec. 

Surveys at 90° Lines of Sight 

Extensive radial surveys were carried out at four axial 
positions in the exhaust plume; 4.6, 10.9, 15.9 and 30.6 cm 
from the anode exit plane, with the line of sight perpendicu- 
lar to the discharge axis. Mirror was adjusted to raise 
or lower the line of sight by intervals of 0.4 cm at the 

plane of the plume axis. At each vertical position the scan 
o 

of the 4880 A All spectral line was recorded from four to 
twelve times, each time during a separate discharge. The 
number of line scans recorded at a given position was dic- 
tated by the reproducibility of the observed line profile. 

In general, higher reproducibility was obtained at the loca- 
tions closer to the anode and near the centerline. 

A typical result, obtained with the line of sight level 
with the axis, 10.9 cm from the anode plane is shown in Fig. 

4-6. The upper trace is the discharge current signature and 
the lower trace is the scan of the 4880 A line of All by the 
Fabry-Perot interferometer. (Two scans of the line during the 
second half of the current pulse are displayed, as discussed 
in section 4-2.) The spectral line exhibits the characteris- 
tic Doppler split resulting from the opposed radial velocity 
components of the two argon jets intersected by the line of 
sight. For clarity the laser reference line was not super- 
imposed simultaneously, but was recorded separately immediately 
prior to the discharge. (It is the small peak midway between 
the two components from the discharge.) It was ascertained 
that the interferometer did not drift in the time between the 
two records. It can be seen that the two Doppler-shifted line 
components from the argon jets are symmetrical about the laser 
line, attesting to the equality of the radial components of 
velocity in the two jets. The second small laser line at the 
position of the next scan of the discharge line at the right 
of the oscillogram was recorded to serve as a relative cali- 
bration of the free spectral range of the interferometer. 
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THE 4880 A All LINE AT 10.9 cm FROM ANODE 


FIGURE 4-6 
AP25 P 395 
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(An absolute calibration of the free spectral range was car- 
ried out with the help of the Steinheil prism spectrograph 
and a high pressure mercury light source, according to a 
technique described by Born and Wolf ^ ^ . The free spectral 

range was set to 0„ 78 A for the bulk of the experiments 
described here. It could not be made smaller because of the 
constraints imposed by the 40 nsec duration of the laser 
pulse. ) 

Figure 4-7 shows a schematic of the exhaust plume with 
the location of the 90 lines of sight used in the survey and 
typical resulting line profiles. Note how the Doppler split in- 
creases with distance from the chamber exit plane, indicating ac- 
celerating radial flow in the argon jets. The blue-shifted compo- 
nents corresponding to the near jet appear always to have a smaller 
intensity than the red-shifted component. The cause for this 
was later found to be due to a slight azimuthal misalignment 
of the MPD accelerator in the vacuum tank, so that the two 
"central" argon jets were actually somewhat tilted out of the 
horizontal plane. The error in measured velocities resulting 
from this is insignificant. The radial velocity components 
of the argon jets are shown plotted against axial distance in 
Fig. 4-8. Each error bar contains the range of velocities ob- 
tained from observed Doppler splits of four to twelve spectral 
line scans at each horizontal line of sight. Because the widths 
of the spectral line components are narrow it can be assumed 
that the spread of radial velocity in each jet is small so that 
the curve shown is a good representation of the development of 
the average radial velocity components in each jet. It can be 
seen that the bulk of the acceleration occurs within about one 
and a half anode orifice diameters (15 cm) from the chamber 
exit plane. 

Figure 4-9 displays the shape of the 4880 A spectral line 
with increasing distance above the axis in the vertical plane 
15.9 cm from the chamber exit. It can be seen that the two 
Doppler-split line components come closer together and eventual- 
ly merge as the optical line of sight approaches the outer 
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edges of the exhaust plume. This behaviour is typical of 
that observed at the other three axial stations and indicates 
the expected decrease in the components of radial velocity pro- 
jected along the line of sight. There is also a noticeable 
decrease in line intensity at the larger radii. 

Radial and Axial Intensity Profiles 

If the plasma is optically thin, the intensity I ( A) 
within a small wavelength interval d A of a spectral line is 
proportional to the number of particles within the volume 
sampled by the line of sight which radiate in this wavelength 
interval. ; Thus the total line intensity = I l(A)dA, 

. . . ''O 

which is simply the area under the line profile, is propor- 
tional to the total number of particles undergoing the spec- 
tral transition giving rise to the spectral line in question. 
The areas under the Doppler-split line profiles of Figs. 4-6 
and 4-7 should consequently be proportional to the sum of the 

total number of particles in the sampled volume radiating at 

o 

the nominal wavelength of 4880 A in each of the two argon 
jets intersected by the line of sight, if the jets are op- 
tically thin to this radiation. 

Experiments were done at the lines of sight 4.6 and 10.9 

cm from the anode plane to determine whether the Doppler-split 
o 

4880 A line is indeed optically thin. The technique consisted 
of placing a concave mirror at twice its focal length from the 
discharge axis in the extension of the optical path and measur- 
ing the increase in intensity. Figure 4-10 shows a schematic 
of this arrangement. If no losses occur at the mirrors or the 
tank windows the observed intensity with the concave mirror 
in place should be double that without the mirror. In actual 
fact there are absorption and reflection losses at the mirrors 
and windows, and these must be corrected for. These losses 
were measured experimentally with the argon laser and a silicon 
photodiode and found to amount to 36%. The MPD arc was fired 
several times, alternately with and without the concave mirror 
in place. The observed sets of line intensities were averaged 
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and compared. The experimental results agreed with the pre- 
dicted increase in intensity within a few percent, well with- 
in the experimental error, thus confirming that the 4880 A 
line is optically thin. 

To obtain radial relative intensity profiles the area 
under each spectral line profile obtained in the radial sur- 
veys was measured (with corrections for background continuum) 
and the results at each distinct line of sight were averaged 
and normalized to the maximum value observed. Figure 4-11 
shows the variation of the relative 4880 A line intensity with 
vertical distance above and below the discharge axis at the 
four axial planes investigated. The off-center maxima corres- 
pond to the upper and lower pairs of argon jets visible in Figs. 
3-9a and 3-9b. The central peak corresponds to the middle pair 
of argon jets. A plot of the loci of the off-center maxima, 
suitably corrected for the 30° azimuthal inclinations of the 
upper and lower argon jets with respect to the vertical axial 
plane, is shown in Fig. 4-12. In effect, these loci, which turn 
out to be straight lines inclined at 15° to the discharge axis 
with a virtual origin at the cathode tip, are reasonable approxi- 
mations to the average flow direction in each argon jet. The 
. , o 

axial variation of the relative 4880 A line intensity in an in- 
dividual jet is displayed in Fig. 4-13. The dependence of in- 
tensity on distance from the anode plane appears to follow an 
inverse square law. Because the jets appear to expand laterally 
in a linear fashion the integrated path length through the jets 
increases linearly with increasing axial distance. Thus, if the 
radiating particle density were to vary as the inverse square of 
the distance from the chamber orifice, such as in source flow- 
like free expansions, one would expect an integrated intensity 
which decreases only inversely with distance. That the intensity 
decreases faster is an indication that the radiating particle 
density falls off more steeply than the inverse square of distance. 

Since the radiating particle density is a strong function of elec- 

(25 2 ^ J 

tron temperature, - 'the above results signify that the electron 
temperature decreases with increasing distance from the chamber 
exit. Further discussion of the electron temperature is given 
in section 4-4. 
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DEVELOPMENT OF 4880 A LINE INTENSITY IN AE JET 

FIGURE 4- 13 
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Survey at 75° Lines of Sight 

A series of Doppler-shift velocity measurements was 
done at five axial positions with the line of sight at an 
angle of 75° to the discharge axis, passing through the two 
diametrically opposed argon jets in the horizontal plane, as 
shown in Fig. 4-14. Assuming as a first approximation, from 
the results of the previous section, that the average jet flow 
angle is 15° with respect to the axis, these lines of sight 
are perpendicular to the far jet and intersect the near jet 
at an angle of 60°. The observed line profile should there- 
fore have a shifted and an unshifted component. This is al- 
most the case, as can be seen in Fig. 4-15, obtained at a 
line of sight which crosses the axis 16.8 cm from the chamber 
exit plane (Line of sight " c" in Fig. 4-14) . The upper trace 
in each oscillogram corresponds to the discharge current his- 
tory, while the lower traces exhibit line profiles scanned 
by the Fabry-Perot interferometer. In each case two successive 
scans of the line were recorded; the comments here apply to 
the first scans on the left. Figure 4-15a shows the relative 
Doppler shift between the near and far jets. The "unshif ted" 
component on the left corresponds to the far jet, intersected 
at — 90°, while the blue-shifted component to its right cor- 
responds to the near jet intersected at ~60 o . Figure 4-15b 
displays the absolute Doppler shift between the line component 
from the near jet and the simultaneously recorded 4879.9 A 
reference laser line. (The laser intensity was made as large 
as possible relative to the discharge intensity to minimize 

distortion of the laser line by the "unshifted" component of 
o 

the 4880 A line from the far jet. This is the reason for the 
different intensity scale in Fig. 4-15b) . The shift between 
the near jet line component and the laser reference is some- 
what less than that between the line components of the near and 
far jets, indicating that the true mass-averaged velocity vec- 
tor of the argon in the jets makes an angle somewhat greater 
than 15° with the discharge axis. 
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In Fig. 4-16 are shown typical absolute Doppler shifts 
between the near jet line component and the laser reference 
line, obtained at the four other 75° lines of sight. The 
vertical intensity scale in each case is different but this 
does not affect the wavelength calibration. At least ten 
such profiles were recorded at each line of sight and the 
respective Doppler shifts averaged. The line profiles in 
Fig. 4-16 appear broader than the ones obtained at the 90° 
lines of sight (Fig. 4-7) because of the failure of one of 
the fine-adjust piezoelectric transducers on the stationary 
interferometer plate mounting. This caused considerable 
difficulty in the adjustment of plate parallelism in the 
interferometer, and the overall instrumental finesse de- 
creased from 25 to about 12 as a result. Fortunately, the 
loss of resolution is not serious. 

Calculation of Jet Velocity and Flow Angle 

Once the velocity components at 75° and at 90° to the 
discharge axis are known at various positions along the argon 
jets, it is relatively straightforward to deduce both the mag- 
nitude and the direction of the average velocity vectors in 
the argon jets as functions of axial distance from the chamber 
exit plane. An initial estimate of jet flow angle and locus 
is required in order to obtain the approximate positions along 
the near jet at which the 75° lines of sight intersect. It 
has already been shown that the peak luminous intensity in each 
jet follows a straight line inclined at 15° to the axis with 
a virtual origin at the tip of the cathode (Fig. 4-12). This 
locus will be used initially to define the position of the 
argon jets. It will turn out that this approximation is good 
enough that no further refinements will be required for the 
calculations outlined below. 

The calculation procedure is as follows: choose, for ex- 
ample, the 75° line of sight labeled "c" in Fig. 4-14, which 
intersects the discharge axis 16.8 cm from the chamber exit 
plane. From the figure it can be seen that it also intersects 
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the near- jet locus at a point 18.5 cm from the chamber exit 
and 6.2 cm from the axis. Figure 4—8 reveals that the average 
radial component of velocity V_ is 5.0 xlO 3 m/sec at this 

IX 

point. The measured Doppler shift along the chosen 75° line 
of sight yields an average velocity V_ = 8.4x10 m/sec. The 
average velocity Vg along the jet and the average flow angle 
oc can be obtained from the following simple relations: 


Solving for o< : 


Thus 


Vi = V s cos ( 75’ - et) 
V R - V s cos ( 30° - *.) 


<x - 


cot r 3.86 ( 

L \ V R 


- 0.965 


’)] 




Vl 

slncx. 


C05 


% 

(75°-oC) 


A simple graphical construction, as shown in Fig. 4-17, can 

also be used. For the example chosen the results are 
4 o 

Vg = 1.47 x 10 m/sec, <X = 20 . Because of the relatively flat 
radial velocity profile (Fig. 4-8) the discrepancy between the 
assumed and calculated flow angles does not warrant recalcula- 
tion of the velocity Vg with a better approximation of c* . The 
errors in V g arising from the original assumption of cx; are 
well within the experimental uncertainties. 

The' results of the velocity and flow angle calculations 
are displayed in Figs. 4-18 and 4-19. Figure 4-18 shows a 
plot of jet velocity V g versus axial distance, x , from the 
chamber exit. The circular data points were obtained from 
Vg = V Q /cos (75°- o<) while the triangular data points were calcu- 
lated from Vg = V R / sin <X , where V R was obtained from the data 
at the four lines of sight used in the perpendicular surveys. 
(Fig. 4-8) The vertical error bars represent the degree of 
uncertainty arising from the spread of the experimental data, 
while the horizontal error bars encompass the probable error 
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in axial position resulting from the original assumption of 
15 flow angle. The average velocity along each jet is seen to 
increase from ^ 8300 m/sec 4.6 cm from the chamber exit to a 
plateau of ~16,500 m/sec 40 cm downstream. The bulk of the 
acceleration occurs over a distance of the order of an anode 
orifice diameter. It is interesting to note that the observed 
final velocity is nearly twice the so-called Alfven critical 
velocity for argon (8700 m/sec), (^,30) which has been recently 
claimed by Malliaris et al. ^ 20 ^ The present result indicates 
that specific impulses as high as about 1600 sec may be at- 
tainable using argon as a propellant, contrary to the limiting 
850 sec suggested by Malliaris A discussion of possible 

acceleration mechanisms which may be responsible for the high 
argon velocities observed here is presented in Chapter 5. 

Figure 4-19 shows the variation of average jet flow angle, 
o< / with distance from the chamber exit. o< increases from about 
17 at x = 6 cm to about 20 at x = 18 cm and remains constant 
at about 20° thereafter. This is a somewhat different result 
from that anticipated from the locus of maximum jet intensity 
and is a result of line of sight integration effects. In 
actual fact, the 15 locus is an accurate representation of 
the local flow angle rather than the average flow angle which 
is weighted to larger values by the higher flow angles present 
at the larger radii 0 

4-4 TEMPERATURE MEASUREMENTS 

In order to determine whether the high flow velocities 
observed in the argon jets are a result of electrothermal con- 
version via the downstream expansion permitted by the MPD ac- 

(13) 

celerator geometry, it is necessary to measure the axial 

profiles of ion and electron temperature in the argon jets. 
Kinetic atom and ion temperatures in plasmas of 1-100 eV are 
generally determinable from measurements of half-widths of 
Doppler-broadened spectral lines, the technique that is used 
here. Electron temperatures can be measured spectroscopically 
by the ratio of line intensities, but only in plasmas which 
are either in full or local thermodynamic equilibrium.^ 25,40 ^ 
Such conditions do not in general exist in the plasma flows 
produced by MPD arcs (2) and electron temperatures thus are 



89 


(2 8 281 

usually measured by means of electrostatic probes. ’ ' 

This section discusses in detail the measurement of argon ion 

temperatures with the scanning Fabry-Perot spectrometer, and 

presents electron temperature data obtained by other investi- 

( 8 281 

gators on the same accelerator. ' ' 


Doppler Width Measurements of Ion Temperature 


The motion of a radiating particle toward or away from 
an observer leads to a Doppler shift in the wavelength of the 
emitted line. In a plasma the random thermal motions of the 
radiating particles cause a Doppler broadening of the lines 
as the overall result. If the velocity component of a par- 
ticle parallel to the direction of observation is v , then 

It 

the wavelength shift is 

a\ ~ ± -q- 


where \ Q is the unshifted wavelength and c is the speed 
of light. If the motion of the plasma particles is purely 
thermal the velocity distribution of the emitters is Max- 
wellian. Thus, the fraction of particles moving in the line 

of sight with velocity components between v n and v + dv is 
(47) p P P 

given by 



where T is 
mass and k 
Vp according 


the temperature of the particles, m is the atomic 
is Boltzmann's constant. By substituting AX for 
to the Doppler shift relation one obtains 


dn 

71 



If the lines are optically thin, the intensities I are propor- 
tional to the concentration of radiating particles n; and, in 
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particular, the intensity emmitted in the interval d( aX ). 
namely I ( aA )d( A A ) is proportional to the fraction of par- 
ticles dn. Thus, for a purely Doppler-broadened line one 
obtains a Gaussian shape for the intensity distribution : ^ 4< ~* ^ 



where 1^ is the total line intensity and AX is the devia- 
tion from line center A 0 - The intensity reaches one-half 
its maximum value when the exponent has the value 1/2. The 
full half-width, which is the width between the two half- 
maximum intensity points follows as: 



(- 


'm c z J 


o 


= 716 x to 7 



where a \ d is obtained in A if X is in A and T is in °K, 
and M is the atomic weight. It is evident that thermal Doppler 
broadening is most pronounced for lines of light elements at 
high temperatures . In the case of the 4880 A line of ionized 
argon temperatures in the range of 10 4 to 10 5 °K result in half- 
widths of 0.06 to 0.19 A. 

Doppler broadening is not the only important line broaden- 
ing mechanism found in plasmas. A variety of other mechanisms 
such as Stark broadening, Zeeman splitting, and non— thermal 
motions such as gross mass motion and microturbulence will 

broaden spectral lines to a degree which depends on the ore- 
,... (25,40,48,49) V 

vailing conditions.' Other effects which act to 

distort line profiles are self-absorption, if the spectral line 
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is not optically thin, and instrumental broadening due 

to the finite resolution of the recording spectrometer. 

It is therefore essential to ascertain experimentally or 
theoretically that either these broadening effects are in- 
significant or that they can be accurately accounted for 
in the data reduction. The importance of each of the above 
mentioned broadening mechanisms will now be examined. 

Stark Broadening? 

Stark broadening of lines emitted by atoms or ions in 
a plasma is a result of the interaction of the radiating par- 
ticles with electric fields originating in the surrounding 

ions and electrons, and is a direct function of the electron 
(25 40) 

density. ' The quantum mechanical calculations which 

are used to determine the relation between the broadening and 
the electron density for a variety of atomic species are dis- 
cussed in detail in Refs. 25 and 40. Recent experimental 

measurements of the Stark broadening of several lines of 
( 50) o 

ionized argon show that for the 4880 A line the Stark 

o 

half-width amounts to ^0,05 A at an electron density of 
16 3 

10 cm . Measurements of the electron density in the ex- 
haust of the MPD arc reveal densities of the order of 10 15 cm 

and less in the regions of the plume downstream of about 5 cm 

(15) 

from the chamber exit. Since the Stark width is directly 

proportional to the electron density for heavy atoms and 
ions (25,40, 50) j_ s eX p ec ted that the Stark width of the 
4880 8 line of All will be less than 0.005 A in those regions 
of the argon jets of interest here. This width is more than 
an order of magnitude less than the estimated thermal Doppler 
widths and may thus be neglected with little resulting error. 
Zeeman Splitting: 

The presence of strong magnetic fields will cause a 
Zeeman splitting of spectral lines and thus alter the line 
profiles. The order of magnitude of the splitting of a line 
of argon of wavelength A (A) is given by^^ 
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where B is the magnetic field strength in Weber/m 2 . The 
highest field strengths found in the MPD accelerator occur 
in the vicinity of the cathode and have a magnitude of about 
0.3 Weber/m (see Chapter 5, Fig. 5-1) leading to a Zeeman 
s Plit of the 4880 A line of — 0.007 A in those locations. 
However, 2 cm downstream of the chamber exit B is down to 
0.01 Weber/m and the Zeeman splitting is thus only about 
2 x 10 A, a negligible amount. Further downstream it is 
even less. 

Microturbulence : 

There is abundant evidence that strong ion turbulence is 
often present in magnetically confined plasmas . * 25 ' 48 ' 49 ^ Al- 
though the exhaust flow of the MPD arc is not in the strict 
sense magnetically confined (except in the vicinity of the 
cathode tip) there is evidence of turbulence. For example 
Langmuir probes and double electrostatic probes display ran- 
dom fluctuations of ion saturation currents which are usually 
ascribed to turbulent fluctuations convected past the probes . 
These may, of course, be manifestations of wave propagation 
offsets but this is not yet clear. Turbulence can be so fine 
grained that the velocity distribution is closely Maxwellian 
on a macroscopic scale while quite large deviations could 
occur on a microscopic scale. If the turbulent motion is 
Maxwellian on the macroscopic scale its contribution to line 
width is determined by^ 25 ^ 




Al 


where v r describes the relative macroscopic velocity com- 
ponent along the line of sight of the various volume elements 
of the emitting plasma and }\ 0 and c have their usual defi- 
nitions. Unfortunately there are no data available on average 
random turbulent velocities in the exhaust of the MPD arc, and 
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until quantitative measurements can be made it will be 
assumed that microturbulent line broading is not signifi- 
cant. * 

Gross Mass Motion: 

Gross mass motion will contribute to line broadening 
if there is a distribution in the velocities of macroscopic 
volume elements of the plasma along the line of sight. 

Clearly, this effect occurs within the argon jets in the 
MPD exhaust, inasmuch as there is a spread in the radial 
velocity components across the jets. The shape of the line 
broadening function arising from this effect depends on the 
radial profile of velocity and local line intensity across 
each jet. These are difficult to evaluate accurately be- 
cause the azimuthal non-uniformity of the exhaust does not 
permit reduction of the spectroscopic measurements to purely 
local values. For present purposes it will be assumed that 
to a first approximation the flow velocity is constant across 

a jet and equal to the average measured velocity, V . Further- 

o 13 

more, it will be assumed that the 4880 A local intensity pro- 
file across each jet has a Gaussian shape. This assumption 
is reasonable because the integrated intensity distributions 
across the upper and lower jets (Fig. 4-11) appear approximately 
Gaussian. This being the case, the broadening function due to 
the distribution of radial velocities in each argon jet will 
also be Gaussian. It remains to determine the radial angular 
spread of the jets at their half-intensity points. From the 

o 

photographs of Figs. 3-9 and 3-10 and the intensity profiles 
of Fig. 4-11 the jets are estimated to have a 10° radial 
divergence angle about the 15° jet locus. With this informa- 
tion and the measured average jet velocities, V g , the com- 
ponents of radial velocity at the inner and outer half-intensity 
points of the jets can be estimated. The difference between 
* 

One method of determining the importance of microturbulence 
is to measure the widths of lines from radiating species of dif- 
ferent charge to mass ratio.( 48 ) If turbulence is significant a 
different temperature will be obtained with each specie. Un- 
fortunately, due to the species separations present in the MPD 
exhaust and the opaqueness of the Fabry-Perot plates to ultra- 
violet radiation from AIII ,it has not been possible to carry out 
such an experiment. 
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these two values at a given axial position gives rise to a 

Doppler shift which provides a reasonable estimate of the 

line broadening due to the gross mass motion. For example 

at x = 30 cm such a calculation yields a broadening of 
o 

~0.04 A which is of the same order of magnitude as the ex- 
pected thermal Doppler widths . 

Line Reversal: 

Distortions to the line profile due to finite optical 

depth of the plasma will also result in a broadening of the 

spectral line. In the present case however, this effect is 

non-existent, since the plasma has been verified to be optically 
o 

thin at 4880 A, (Section 4-3) 

Instrumental Broadening: 

The Fabry-Perot interferometer, like all other spectrometric 

instruments, has a finite resolving power and consequently 

(42 43 45) 

leads to a further broadening of the spectral lines. ' ' ' 

The broadening mechanisms in the Fabry-Perot spectrometer are 
well known (Appendix B) and in general result in an instru- 
mental function which very closely approximates a Gaussian 
shape. This is fortuitous as it simplifies the unfolding of 
the instrumental function from the Gaussian Doppler profile. 

The folding of two Gaussian shapes results in a Gaussian pro- 
file with half-width equal to the square root of the sum of 
the squares of the individual half-widths A 5 "*"' The half-width 
of the Fabry-Perot can be easily calculated from its total 
finesse F T and its free spectral range aAj . As used in the 

experiments described here, the interferometer has a finesse 

, o 

of 25 and a free spectral range of 0.78 A. Therefore the in- 
strumental half-width is given by aA.^ = A\ f/F T = 0.031 A. 

Experimental Results 
o 

The 4880 A spectral line profiles of All obtained in the 
90° line of sight survey were used to calculate the argon ion 
temperature from line widths. Only the data obtained at the 
lines of sight level with the axis, intersecting the two argon 
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jets in the horizontal plane, were used in order to obtain 
maximum separation of the two Doppler shifted line components. 
The estimates of other line broadening effects discussed in 
the previous sub-section have shown that Stark, Zeeman and line 
reversal effects can be neglected compared to expected Doppler 
widths. Microturbulence broadening could not be evaluated be- 
cause of lack of pertinent data. However, the broadening 
caused by the spread of radial velocities across each jet and 
the instrumental broadening were shown to be significant and 
thus were taken into account in the reduction of line widths . 


The observed profile of each line component is a convo- 
lution of the Doppler profile of half-width aA d , the veloc- 
ity spread function of half-width aA-_# and the instrumental 


v' 


function of half-width aA^* Since the latter two are here 

assumed to be Gaussian and the Doppler profile is Gaussian, 

(51 52 ) 

the total measured half-width is simply: ' ' 


a A T = (a A / +■ 




Thus the true Doppler width is: 


A D — ( a\ j- — aX v — AXi ) 


'/2 


and the argon ion temperature is given by 

,13 f A\ 


T = 7.8 x 10' 


W 


= 4880 A becomes 


K 


T = 3.28 x 10 6 A A d 2 °K 

= 2.83 x 10 2 AA^ 2 eV 


where AA d is in A. It must of course be emphasized that 
due to the line of sight integration effects the temperature 
so calculated is only an average across the width of the jet. 
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The results are shown in Fig. 4-20 with the ion tempera- 
ture T^ in the argon jets plotted against axial distance from 
the chamber exit plane. The circled points are the averages 
obtained from all the line widths measured at each line of 
sight and the error bars encompass the spread in measured 
total widths and the uncertainty in the radial velocity dis- 
tribution at each line of sight. It should be mentioned here 
that even under ideal conditions it is difficult to obtain an 
accuracy in the measured temperatures that is substantially 
better than 20% because the line widths themselves are dif- 
ficult to read to an accuracy better than 10%. In the present 
case .the lack of knowledge of the true radial velocity distri- 
bution function and the neglect of microturbulence effects in- 
creases the uncertainty. In general, however, since most 
processes tend to increase the line width rather than decrease 
it, the results shown in Fig. 4—20 may be considered at worst 
to be upper limits of the ion temperature. 

The best line drawn through the circled points indicates 
that the ion temperature follows an inverse power law, 

-3/4 

T i~ x * decreasing from ~ 1.8 eV at x = 4.6 cm to ~ 0.5 eV 

at x = 30.6 cm. Extropolation back towards the chamber pre- 
dicts values of 5—6 eV in the vicinity of the chamber exit. 

The observed ion temperatures are significantly lower than had 
been estimated previously on the basis of spectrographic ^ 18 ^ 
and Langmuir probe data. ( The overestimate in the previous 

spectrographic experiments v ' is now known to have been a re- 
sult of neglecting the Doppler split in the spectral lines due 
to the discrete argon jet structure in the exhaust plume. (It 
is interesting to note that if it is assumed that the exhaust 
flow behaves as a free isentropic expansion so that /D ~ x -2 

i » 

and T ~r , where jO is the mass density and Y is the ratio 

of specific heats, the fact that T.^x -3 / 4 predicts ^«1.38 

which agrees with values deduced by other means in similar en- 

(13) 

vironments, ) The decrease in ion temperature down the plume 
simultaneously with the increase in flow velocity suggests that 
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directed kinetic energy is being recovered from random 
thermal energy. Whether the change in ion temperature is 
sufficient to explain the observed increase in velocity will 
be examined in Chapter 5. 

Electron Temperatures 

Very few data on electron temperature profiles in the 
exhaust plume exist at this time. Turchi^'^ has measured 
electron temperatures with Langmuir probes at several lo- 
cations inside the arc chamber and just outside the exit 
plane. His results indicate that the electron temperature 
is essentially constant throughout the regions investigated 
and has a value between 1.5 and 2 eV. Using double electro— 
static probes Boyle v has recently measured the electron 
temperature on the locus of maximum luminous intensity in 
one of the argon jets, 25.4 cm from the chamber exit, and 
found a value of about 0.8 eV. Comprehensive radial and 
axial surveys have not yet been made, but indications are 
that the electron temperature on the axis, at least at 
x = 25.4 cm, is about the same as in the argon jets, and that 
it xncreases axxally towards the exit plane. ' It would 
thus appear that the electron temperature is radially quite 
uniform everywhere and that it decreases slowly with in- 
creasing distance from the chamber. 

4-5 IMPURITY VELOCITIES AND TEMPERATURES 

Unfortunately, a spectroscopic measurement of the veloc- 
ities and temperatures of the products of ablation (CII, Nil, 
Oil, HI) in the core of the exhaust and in the gaps between 
the argon jets has not proved feasible with the Fabry-Perot 
system. This is because of the lack of suitable isolated 
spectral lines of these impurities within the useful spectral 
range of the reflective coatings on the Fabry-Perot plates. 

The impurity lines tend to crowd together at separations of the 
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order of the free spectral range or less and would thus be 
impossible to identify in the Fabry-Perot scan. Attempts 
at measuring impurity velocities have been made by using 
the Steinheil spectrograph, v but the dispersion of this 
instrument is not high enough even in the blue end of the 
spectrum to yield accurate quantitative results. However, 
the results of that investigation do indicate that the 
velocities in the core flow are considerably higher than 
in the argon jets, possibly as high as 4 - 5 x 10 4 m/sec. 
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CHAPTER 5 / 

ARGON ACCELERATION MECHANISMS 

5-1 INTRODUCTION / 

The experimental velocity measurements described in 

the preceding chapter indicate that almost half the final 

argon exhaust velocity of 1.65 xlO 4 m^sec is achieved in 

an acceleration region downstream of /the arc chamber exit. 

It has variously been suggested that /acceleration in the 

exhaust downstream of the chamber maj be a result of 

fringing electromagnetic fields lor of recovery of 

directed kinetic energy from thermal energy via gasdy- 
• . (13) 

namic flow expansion. A momentum exchange between the 

core flow of ablation products and the argon jets is also 
a possibility. 

The flow dynamics of the MPD exhaust plume would be 
very complicated even without the complex species structure 
the exhaust is now known to have. Highly non-equilibrium 
conditions characterize the entire flowfield, and species 
in various stages of ionization are present. Close to the 
chamber exit are found regions of significant current con- 
duction and self-magnetic fields. The exhaust flow, once 
outside the arc chamber, is unbounded and far enough down- 
stream is expected to exhibit behaviour typical of free jet 
expansions. Theoretical investigations of the flow of high 
enthalpy plasmas into a vacuum have been carried out by some 
investigators, but usually with many simplifying assumptions . ^ 54 ^ 
In the present case there is the added complication of the 
radial and azimuthal species inhomogeneity and resulting non- 
uniformity in the atomic weight of the particles in the flow- 
field. A complete theoretical analysis of the MPD exhaust 
flow is not within the context of this dissertation. In- 
stead simple arguments are developed for the evaluation of 
the relative roles of the three possible argon acceleration 
mechanisms pointed out above. 
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5-2 FRINGING ELECTROMAGNETIC FIELDS 


An estimate of \the role of the j x B body forces in 
the flow acceleration downstream of the anode plane can 
be obtained from a simple comparison of the magnetic, gas- 
kinetic, and dynamic pressures in the acceleration region, 
which lies approximately between x as 5 cm and x w 15 cm. 

Magnetic Pressure \ 

(4) \ 

Clark' has obtained maps of the magnetic field strength 
in and around the discharge chamber of the MPD accelerator. 
Figure 5-1 displays the contours of self-induced magnetic 
fields obtained at a curr'ent of 16 kA and an argon mass 
flow of 6 g/sec. The actual field is, of course, azimuthal. 
Superimposed on these contours is the approximate location 
of an argon jet. It can be seen that at x » 5 cm in the jet 
B « 3 x 10 Weber/m , which is about two orders of magni- 
tude less than in the vicinity of the cathode. The magnetic 
pressure at this axial location in the jet is thus^ 14 ^ 


= e_ * 4 

Zfo * 


Ntn/n/ 


where y x q = ATT x 10 henry/m is the magnetic permeability 
of free space. Downstream of this location B decays rapidly 
and thus so does p^ — at x = 6 cm p is nearly an order of 
magnitude smaller than at x = 5 cm. 


Gaskinetic Pressure 


To estimate the gaskinetic (static) pressure it will 
be assumed that the electron and ion temperatures at x = 5 cm 
are equal and have a value of ~ 1.7 eV (2 xlO 4 °K) as 
measured for the argon ions. (Chapter 4, Sect. 4-4). 
Particle densities are estimated on the basis of electron 
density measurements obtained from Stark width measurements 
of the Hp line. ^ 15 ^ Figure 5-2 reproduces a map of electron 
density contours obtained at 16 kA and 6 g/sec. The approx- 
imate location of an argon jet is superimposed. At x = 5 cm 
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the electron density n has an average value of about 
15 -3 . e 

2 x 10 cm in the argon jet. Although nonequilibrium con- 
ditions generally prevail throughout the flow-field, in 
order to approximate the heavy particle number density it will 
be assumed that Saha equilibrium exists at this location in 
the jet. At the conditions of T * 2 x 10 4 °K and n 2xl0 15 crrf 3 , 

/ 55 \ ® 

equilibrium thermodynamic tables v ' indicate that the den- 
sity, n^ , of singly ionized argon is ~ 7 x 10^ cm -3 and the 
density, n^ , of doubly ionized argon is ~ 8 x 10 44 cm -3 . The 
number densities of neutral atoms and triply ionized atoms 
are negligible. 

The static pressure, p , is calculated from the ideal 

s 

gas law: 

p s = * n xl + nz)kT 

3 2 

» 10 Ntn/m 

where k is Boltzmann's constant, 1.38xl0 -23 joules/°K. 

Dynamic Pressure 

Although experimentally measured radial profiles of 
impact pressure at several axial stations are available/ 1 ^ 
there are no data at x = 5 cm. The dynamic pressure can be 
obtained by interpolating between existing profiles and 
subtracting the calculated static pressure; however, both 
the axial and radial impact pressure profiles are steep 
near x = 5 cm making interpolation difficult and inaccurate. 

A simple way to calculate the dynamic pressure, p, , is through 

. 2 ° 
the relation p^ = pu , where jO = rru (n^ + n^) is the mass 

density, m^ is the atomic mass of argon, and u is the flow 

velocity. Using the measured average argon velocity - 8.5 xlO 3 

m/sec - at x = 5 cm and the previously calculated particle 

densities at this location one obtains 

frjL ** ^ Ntn/in 
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Comparison of Pressure Components 

The preceding calculations show that the magnetic 
pressure is more than two orders of magnitude less than the 
gaskinetic pressure and three orders of magnitude less than 
the dynamic pressure at x = 5 cm, i.e: 

<< << Af 

It can be shown that with increasing x, p m decreases faster 

than p or p, , so that this relation remains valid for x>5 cm. 
d 

The smallness of the magnetic pressure term in comparison 
with the other pressure components clearly indicates that 
electromagnetic body forces play a negligible role in the 
acceleration of the argon plasma in the exhaust in the region 
of interest. 

5-3 ELECTROTHERMAL CONVERSION 

It is conceivable that the observed acceleration in the 
argon jets is a result of the conversion of thermal energy to 
directed kinetic energy via the free flow expansion processes. 
In the preceding chapter it was ascertained that both the ion 
and electron temperatures decrease with increasing distance 
from the arc chamber exit. The ion temperature, T i , was shown 
to have a functional dependence on distance, x , of the form 
T . ~ x - . The exact functional behaviour of the electron 

i 

temperature, T g , is not known at this time. Therefore a simple 
approximation will be used: the known value of T e « 0.8 eV at 
x = 25.4 cm and the assumed value of T e « 1.7 eV (based on the 
results of Turchi^) at x = 5 cm are plotted on full logarith- 
mic graph paper and the function is approximated by a straight 
line between the two points. The slope of this line turns out 
to be very nearly -1/2, and thus the approximate dependence of 
electron temperature on distance is T q ~ x 2 , with T g « T^ at 
x = 5 cm. 

A simple stream tube model can be used to estimate whether 

the decrease in T. and T beyond x = 5 cm is sufficient to 

l e 

account for the increase in flow velocity along the argon jets. 
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It is assumed that the flow in each argon jet is one-dimen- 
sional, i.e: that each argon jet is a stream tube with uniform 
average flow properties at any cross section. Interactions 
with the adjacent flow of ablation products are neglected as 
are radiation and conduction effects. In addition, it is 
assumed, as was done previously, that the flow is in Saha 
equilibrium at x = 5 cm. No assumptions about its state 
downstream of this location need be made at this time. Con- 
sidering the true complex nature of the flow-field the fore- 
going assumptions are admittedly crude; however, enough of 
the essential physics is retained that the results will be 
meaningful in an approximate sense. 

The one-dimensional energy equation for continuum 

adiabatic flow in an argon jet, under the above assumptions 

• (56) 

is 

2 , 

h T = h + j£- - const. 

where h^ is the total enthalpy, h is the static enthalpy, and 
u is the flow velocity at a given cross section of the jet. 

Energy Inventory 

The static enthalpy of a gas is given by h = e + p/p 
where e is the internal energy, p is the static pressure, 
and p is the mass density. The internal energy consists of 
the random thrmal energy of translation, e , of the ions and 
electrons, the enrgy of excitation, e , of the atoms and 
ions, and the energy, e^ , of ionization of the atoms and 
ions. Thus, 

A - e +■ fir + e. + jfL 

€ X 

Using the experimental data under the assumption of Saha 
equilibrium the components of h are easily calculated with 
the help of tables of thermodynamic properties of argon^ 5 ^. 

At x = 5cm the results are as follows, expressed in units of 
electron volts per heavy particle; 
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Translational energy: 

e tr 

~ 6.6 

eV 

Excitation energy: 

e e 

* 0.9 

eV 

Ionization energy: 

e . 
i 

« 31 

eV 

Flow work: 

p/f 

« 4.3 

eV 


Thus the static enthalpy per heavy particle at x = 5 cm is 

h » 43 eV 
3 

At this location u « 8.5x10 m/sec, and the directed kinetic 
2 

energy, ^ ixku , per argon ion is ~ 15 eV. Consequently, the 
total enthalpy per argon ion at x = 5 cm is 

/)_ = h + « 58 eV 

T 2 

4 

The final exhaust velocity of the argon jets is 1.65 x 10 
m/sec, which corresponds to a kinetic energy per heavy par- 
ticle of 57 eV. It is apparent, therefore, that in the ab- 
sence of other acceleration processes the entire thermal 
energy content of the flow in the argon jets present at 
x = 5 cm must be recovered as directed kinetic energy to 
produce the observed final argon exhaust velocity. 

This, however, does not appear to be the case. Dis- 
charge spectra obtained at several axial locations in the 
exhaust plume (Chapter 3, Section 3-3) reveal no evidence 
of neutral argon, whose spectral lines should be visible 
if there is significant recombination of the singly ionized 
argon. All. Furthermore, All radiation persists far down- 
stream. These facts immediately imply that 15.8 eV (the 
first ionization potential of argon) per heavy particle are 
not available for conversion to directed kinetic energy. 

If only this energy of first ionization is non-recoverable 
as streaming energy, and if all the other components of the 
enthalpy are, it turns out that the available enthalpy at 

4 

x - 5 cm would lead to a final argon velocity of ~1.4xl0 
m/sec, which is within 15% of the measured value, a tolerable 
discrepancy in view of the rather drastic assumptions made 
in the analysis so far. If this result is to hold it must 
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first be determined whether the energy of second ionization 
is recoverable. 

Recombination of AIII to All 

Although there is no evidence of All recombination, 
indications are that there may be some recombination of AIII 
to All. For example, the 4880 S All line intensity mea- 
surements in Chapter 4 reveal a line of sight integrated 

_2 

intensity which varies as x along each argon jet. Because 

-t 

of the axial decrease in electron temperature as x 2 and 

the exponential dependence of excitation on this temperature, 

_2 

the line intensity would decrease much faster than x 
unless the upper level of the 4880 X spectral transition 
were populated by cascading transitions from higher energy 
levels of All. Such cascading would be a natural consequence 
of three-body recombination from the next higher level of 
ionization, namely AIII. 

Under the conditions of temperature and particle density 

prevailing in the exhaust plume of the MPD arc the dominant 

recombination mechanism is indeed that of three-body 

. (57-69) 

electron-electron-ion collisions. The energy re- 

leased by the recombination of an ion of AIII to an ion of 
All is simply equal to the second ionization energy of 
argon. Some fraction of this energy is emitted as radiation 
and the rest is transferred to the electron gas via super- 
elastic electron-ion collisions. The relative magnitude of 
the two effects depends largely on whether or not the plasma 
is optically thin to the emitted radiation. ^ 63,68 71) 
Because the treatment of even one-dimensional recombining 
flows with only one ionic specie present is a matter of con- 
siderable complexity, a full analytical solution of the prob- 
lem has not been attempted. Instead, a semi-empirical ap- 
proach based on the available experimental data has been 
used to estimate the fraction of recombination of AIII to 
All whithin the acceleration region in the argon jets. The 
details of the calculation are discussed in Appendix C. 
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The results of the calculations indicate that only 
about 7% of the second ions present at x = 5 cm recombine 
downstream. If the entire energy of this fractional recom- 
bination were recovered as random thermal energy of the 
electrons and ions, and subsequently as streaming energy, 
only about 1 eV per heavy particle would be regained - an 
insignificant amount. Actually, under the conditions preva- 
lent in the argon jets, more than half the recombination 
energy is lost as radiation, as can be ascertained through 

(73) 

calculations similar to those of Fraser et al. 

It is clear, therefore, that neither the energy of first 

or second ionization is recovered in the flow expansion, so 

that frozen flow conditions prevail in the argon jets, and 

very likely in the entire exhaust plume. Only the enthalpy 

of the random thermal motions of the electrons and ions is 

available. The total velocity increment resulting from the 

complete conversion of this enthalpy (^11 eV) into streaming 

3 

motion downstream of x = 5 cm is ‘-'2.7 x 10 m/sec, which is 

3 

less than 35% of the observed increment of ~8xl0 m/sec. 

One factor which has not been taken into account in the 
calculations above is the contamination of the argon jets by 
the ablated impurities. In Chapter 3 it was seen that while 
the heavy impurities such as carbon, oxygen and nitrogen do 
not appear to penetrate the argon jets to a great extent, 
the hydrogen does as a result of its high mobility. The pre- 
sence of hydrogen in the argon jets effectively lowers the 
average atomic weight of the plasma in the jets, allowing 
gasdynamic expansion of the flow to higher final velocities 
than pure argon can achieve. However even if there were a 
hydrogen atom or ion for every argon ion in the argon jets, 
thus lowering the average atomic weight from 40 to about 20, 
gasdynamic expansion would yield a velocity increment of 

3 

5 x 10 m/sec, which is still considerably short of the ob- 

3 

served 8x10 m/sec. Such a concentration of hydrogen m 
the argon jets would require an ablation rate of more than 
one-third the injected mass flow rate of argon. In fact. 
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indications are that the ablation rate does not exceed about 
15% of the input mass flow at the 16 kA, 6 g/sec operating 
condition, implying that the average atomic weight in the 
argon jets is closer to 30 than to 20. Thus the velocity 
increment through electrothermal conversion is at most only 
half the observed value. 

5-4 MOMENTUM TRANSFER FROM THE ABLATED FLOW 

A possible explanation for the large increase in velo- 
city downstream of x = 5 cm in the argon jets can be found 
in the centerline time-of flight velocity profile reported 
in Ref. 11, which is shown here in Fig. 5-3. According to 
the results of Chapter 3 this profile corresponds to the 
velocities of the ablation products along the centerline, 
since there is no argon in the core of the exhaust plume. 

The velocity of the impurities on the centerline is consi- 
derably higher than that of the argon in the jets, a conse- 
quence, most likely, of the much lower average atomic weight 
of the ablated products ( M at)1 S! * 6 -7 vs. M = 40 for argon). 
There is, of course, the possibility that an ion-acoustic 
wave component is superimposed on the convective velocity 
of the ablated flow. Unfortunately, as pointed out in 
Chapter 4, a spectroscopic measurement of the centerline 
velocity profile was not possible; however at this stage 
it is not so much the exact magnitude of the centerline 
velocity profile that is important, but rather the curious 
dip in the profile in the region between about x = 5 cm 
and x = 15 cm. It has been suggested that this may be a 
result of the superposition of a decreasing wave velocity 
profile and an increasing convective velocity profile. 

However, it is interesting to note that this pronounced 
dip occurs over precisely that axial interval over which 
the argon in the jets undergoes the bulk of its downstream 
acceleration. It is tempting to conclude from this that there 
is a transfer of axial momentum from the fast core flow of 
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ablation products to the slower argon jets in the axial inter- 
val where the centerline profile shows the velocity defect. 

The increase in centerline velocity beyond the dip could 
then be interpreted as due to an effective decoupling of the 
ablated and argon jet flows resulting from the rapid drop in 
particle number density caused by the free flow expansion. 

Once decoupled, the ablation products and the argon would 
continue to expand independently, attaining final velocities 
which are functions of their respective atomic weights. 

Indications are that the flow in the exhaust plume is 
turbulent rather than laminar, SQ .^at the momentum 

transfer mechanism may be turbulent shear stress between the 
impurity and argon jet flows. The analysis of the interaction 
of two concentric supersonic turbulent jets of different 
atomic weight is in itself formidable. ' In the present 
case the problem is complicated by the azimuthal non-unifor- 
mity of the argon and impurity flows and the resulting non- 
uniformities in atomic weight and velocity profiles. The 
other complications accompanying freely expanding flows of 
high enthalpy plasmas, as discussed earlier, are also present. 
Furthermore, in the axial interval of interest there are no 
detailed experimental radial or azimuthal profiles of impu- 
rity and argon velocity and concentration available with 
which to compare theoretical predictions. It is pointless, 
therefore, to attempt a detailed analytical treatment at 
this time. 

At present the momentum transfer hypothesis is at best 
speculative because of its reliance on the time-of-f light 
velocity measurements, which may not reflect the true streaming 
velocity of the ablated flow on the centerline. The proposed 
mechanism of argon acceleration must therefore remain purely 
conjectural until the problem of insulator ablation is elimi- 
nated. Spectroscopic velocity measurements on an uncontami- 
nated argon exhaust should resolve whether the presently ob- 
served argon velocities are largely the result of interaction 
with the flow of ablation products, or of some other as yet 
unspecified acceleration mechanism. 
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CHAPTER 6 

SUMMARY AND CONCLUSIONS 


The photographic and spectroscopic investigations of 
the structure of the exhaust plume and of the exhaust 
velocities and temperatures of argon have provided impor- 
tant new information on some of the physical processes 
occurring in the quasi-steady MPD accelerator. A hither- 
to unsuspected azimuthal and radial structure in the 
spatial distribution of the injected argon has been found 
to exist in the discharge region and exhaust plume, under 
the standard operating conditions of 16 kA current and 
6 g/sec argon mass flow. Singly ionized argon, the dom- 
inant radiating specie of the propellant, is distributed 
in six discrete jets which are azimuthally in line with 
the six mass injector ports and retain their identity 
for several anode orifice diameters downstream of the 
chamber exit. The gaps between the jets, which have been 
shown to be devoid of any specie of argon, are occupied 
by impurities - carbon, hydrogen, oxygen - ablated from 
the Plexiglas insulator end wall of the arc chamber. 
Operation of the arc at overfed conditions reduces the 
amount of the ablated material in the exhaust and causes 
the argon distribution to become more uniform. Operation 
at starved conditions substantially increases the erosion 
rate of the Plexiglas insulator and exacerbates the com- 
plex species structure in the exhaust. 

The insulator ablation, which is responsible for the 
complex exhaust structure, occurs over regions of argon 
starvation between the injectors. Kerr-cell photographs 
of the early stages of the discharge indicate that the on- 
set of ablation is a relatively slow process compared to 
current pattern stabilization, requiring about 150 fasec to 
achieve a steady state, as compared to 30 |asec for the 
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current. The species structure of the exhaust reflects 
this growth of ablation, and does not develop its full 
structure until about 150 pisec into the current pulse, 
simple calculations indicate that the replacement of the 
Plexiglas insulator by a refractory material such as quartz 
will not entirely eliminate the erosion problem but will 
reduce it significantly and confine it to the immediate 
vicinity of the cathode base, where the radiative and con- 
ductive heat fluxes to the insulator are very high. 

Extensive time-resolved spectroscopic measurements 
of the velocity and temperature of singly ionized argon 
in the jets have been carried out by means of a versa- 
tile scanning Fabry-Perot spectrometer system. The 
average argon ion velocity in the jets has been found to 
increase monotonicaliy from ^ 8.5 x 10 m/sec 5 cm from the 
chamber exit, to a plateau of ~ 1.65 x 10^ m/sec about 40 cm 
downstream. This final velocity is considerably in excess 
of the Alfve'n critical speed for argon, indicating the 
possibility of attaining much higher specific impulses with 
this propellant than the limits claimed by other investiga- 
tors.^ 30 ^ Measurements of spectral line widths, suitably 
corrected for broadening mechanisms other than the thermal 
Doppler effect, have revealed an argon ion temperature pro- 
file in the jets which downstream of x = 5 cm scales as 
-3/4 

x , where x is the distance from the chamber exit. 
Typical temperatures are 1.7 eV at x = 5 cm and 0.5 eV 
at x = 25 cm. 

An evaluation of possible acceleration mechanisms, 
based on these and other data (such as magnetic field pro- 
files and electron temperatures) has indicated that neither 
body forces from fringing electromagnetic fields nor simple 
electrothermal conversion of the random thermal energy of 
the argon plasma are sufficient to account for the high 
velocities attained by the argon in the jets. Calculations 
of the recombination rate of doubly ionized argon have re- 
vealed that less than 10% of the All I present at x = 5 cm 
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recombines downstream, and the spectroscopic studies have 
shown no evidence of recombination of the singly ionized 
specie. These facts indicate that very little ionization 
energy is recovered as streaming energy. 

A comparison of a centerline time-of-f light velocity 
profile with the velocity profile in an argon jet 
suggests that a transfer of momentum from the fast core 
flow of ablation products to the slower argon jets may 
be the principal argon acceleration mechanism. If this 
is indeed the case it might be expected that in the absence 
of the light impurities the argon exhaust velocity would 
be lower than that observed in the present investigation. 
However, it may well be that some other acceleration mecha- 
nism not yet considered may be responsible for the high 
observed argon exhaust velocity. 

It is clear that the ablation of the insulator end 
wall of the arc chamber must be eliminated or at least 
drastically reduced as a prelude to further research and 
development of the quasi-steady MPD accelerator. It has 
been recently determined in this laboratory' ' that 
this can be done through a suitable choice of injector 
geometry and insulator material. The use of annular in- 
jector slits, rather than discrete orifices, has been 
found not only to significantly reduce the rate of abla- 
tion but also to give rise to an azimuthally uniform dis- 
tribution of the injected argon in the exhaust plume. 

The latter is a fortuitous result, as it will permit the 

reduction of future spectroscopic observations to local 

( 16 ) 

properties via Abel inversion techniques; ' Velocity 
and temperature measurements on an uncontaminated axially 
symmetric exhaust plume should resolve the nature of the 
argon acceleration processes present downstream of the 
discharge chamber. 
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APPENDIX A 
ABLATION PHENOMENA 


A-l INTRODUCTION 

It is evident from the discussions in Chapter 3 that 
the ablation of the Plexiglas end-wall of the arc chamber 
has a profound effect on the structure of the exhaust plume 
of the MPD arc. Time-resolved Kerr-cell photographs of the 
arc chamber interior taken during the first 160 psec of the 
discharge at 16 kA and 6 g/sec, indicate that the onset of 
the insulator ablation is a relatively slow process requiring 
about 150 psec for stabilization, compared to about 30 psec 
for current pattern stabilization. A sequence of such photo- 
graphs obtained at 13,21,40,100, 132 and 150 psec into the cur- 
rent pulse is shown in Fig. A-l. Note the radial and axial 
dimensions of the luminous layer of ablation products at the 
surface of the end-wall. At 13 psec the luminosity at the 
wall is just visible near the cathode base. At 21 psec it is 
still within the radius of the injectors. At 40 psec the 
layer around the base of the cathode extends axially halfway 
to its tip and radially to about the injector orifices. Be- 
tween t = 40 psec and t = 100 psec the luminous layer grows con- 
siderably, developing its characteristic star-shaped pattern. 

The bright scalloped area surrounding the cathode is the "core" 
of the star-shaped layer, while the fainter extensions are the 
arms viewed side-on. (Cf. Fig. 3-12) As time increases, so 
does the extent of this layer of ablation products. At 150 
psec the arms of the star reach the cylindrical side-wall and the 
downstream "boundary' r reaches the level of the cathode tip. The 
l a ysr maintains this visual configuration throughout the remainder 
of the current pulse. This growth of ablation can be explained by 
a simple qualitative physical model. After breakdown and current 
stabilization the self— field body forces deplete argon from 
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the regions between the injectors at the Plexiglas wall, 
forcing the current pattern near the surface to ablate the 
insulator material in order to satisfy the local mass flow 
requirements. Near the Plexiglas surface the current dens- 
ity is essentially radial, and hence inversely proportional 
to the radius. it is thus most intense at the cathode base. 
Consequently, the initial argon fill in the chamber is most 
quickly depleted in this region and ablation of the end-wall 
first begins at the cathode base. The radius of argon de- 
pletion increases with time, eventually reaching the side 
walls of the arc chamber. At this time, about 150 psec in- 
to the discharge, the total ablation rate reaches a steady 
state. Because the argon enters the chamber through injec- 
tors extending 6 mm downstream of the Plexiglas insulator, 
very little of it reaches the surface between the injectors 
and consequently the discharge remains starved there. At the 
injectors, on the other hand, there is sufficient argon and 
ablation is much less severe. 


A-2 POSSIBLE ABLATION PROCESSES 

Insulator ablation phenomena similar to those reported 
here have been observed in a variety of different electric 
discharge devices such as 3— pinches and electromagnetic 

11 TO TO on\ 


shock tubes. <76.77.78.79,80) 


Various heating mechanisms for 


the insulator walls have been proposed. Keck ^ 76 ^ working 
with a magnetic annular shock tube has suggested heavy par- 
ticle heat conduction; Colgate ^ 7 ^ , studying pinch dis- 
charges, has proposed electron heat conduction (taking the 
gas to be fully ionized right up to the wall); Ashby^ 78 ^, 
also working with a pinch device, has felt that radiation 
would be the dominant heat transfer process; finally. Workman 
in an investigation of ablation effects in a magnetic piston 
shock tube has suggested the possibility of electron-ion re- 
combination at the insulator surface as an important heating 
mechanism. It is likely that all the above processes play some 
role in insulator ablation. The dominant mechanism probably 
depends on the nature and geometry of the discharge. 


( 79 ) 
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In the MPD arc under investigation here it is believed 
that conductive heat transfer from a thin layer of resis- 
tively heated plasma adjacent to the insulator surface plays 
a significant role in the ablation process. The current 
density at the surface is apparently much higher than any- 
where downstream at a given radius. For example, Fig. A-2 
shows the axial variation of the radial component of steady- 
state current density at a radius equal to that of the cath- 
ode base. These data were reduced from detailed current con- 

(18) 

tours similar to the ones shown in Fig. 2-6b. . It can 

be seen that there is an abrupt rise in the current density 
to more than 2000 A/cm within a layer about 2 mm thick next 
to the end wall. At such current densities the power dis- 
sipated by ohmic heating of the plasma in the thin layer is 
considerable. It is known from plasma potential and current 


conduction patterns that the current conduction near the 
Plexiglas insulator is largely scalar , therefore the 


plasma resistivity 


from 


(14) 


can be calculated approximately 



G.S3 x 10* 

jVe 


ohm- cm 


where T is the plasma temperature and £nj\. is a weak func- 
tion of density and temperature, equal to about 6 for 

15 -3 

T « 1 eV and n *+ 10 cm . For a temperature of about 

/ Q \ © 2 

1.5 eV v , ^ « 1.5 x 10 ohm-cm. Thus the power dissipated 

at the base of the cathode, where j sw 2000 A/cm 2 , is 

P = ^ j Z ~ 6 X 10* W/cm 2 

Because the current varies inversely with radius the ohmic 
heating is proportional to the inverse square of the radius 
in the thin layer at the wall, as shown in Fig. A-3, and is 
thus down to about 9x10 W/cm at the radius of the injectors. 
Some fraction of the heat generated is convected downstream 
by ^xB body forces, and the rest is conducted into the wall 
raising its temperature to the ablation point. The relative 
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magnitudes of these fractions is not known at present. A com- 
plicating factor is that the intense radiation emanating from 
the small volume of continuum emission just downstream of 
the tip of the cathode appears also to play a significant role 
in the ablation of the Plexiglas. This can be seen in 
Fig. A-4 which shows a photograph of the Plexiglas plate re- 
moved from the chamber but with the cathode and the injectors 
still in place. Note the distinct light "shadows" on the sur- 
face around each injector on the sides farthest from the 
cathode. A small light bulb placed near the tip of the cath- 
ode, which simulates the intense continuum emission region, 
is found to cast shadows of the injectors which exactly fill 
these light areas. Examination of the surface in these 
"shadows" reveals only slight damage as compared to other 
areas of the Plexiglas plate. The injectors themselves show 
an erosion pattern which is also consistent with radiation 
heating from the intense source off the cathode tip. No de- 
tailed quantitative studies of the ablation process in the 
MPD arc have yet been made and consequently it is not possible 
to assess the relative roles of heat conduction and radiation. 
It is interesting to recall however, that there is residual 
ablation even when the discharge is highly overfed, which 
suggests that ablation under such conditions may be due to 
radiative effects exclusively. 

A-3 EFFECTS OF INJECTION GEOMETRY AND INSULATOR MATERIAL 

It appears from the foregoing discussions that 
the mass injection geometry has a strong influence on the 
insulator erosion, and that to minimize ablation the mass flow 
of argon should be matched locally to the current density 
everywhere on the surface of the Plexiglas end plate. Re- 
cent investigations lend support to this contention.^ 28,81 ) 
Various injection geometries are in the process of being in- 
vestigated at the present time. In one, the propellant is 
injected radially from each injector head, to provide a film 
of argon over the entire surface of the insulator; ablation 
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is noticeably reduced. Other injection geometries such as 
annular slits and radial injection from the base of the 
cathode meet with even greater success in reducing abla- 
tion. A more sophisticated injection technique that 

may be useful is a porous end wall. 

The use of other more durable insulator material is 
also being considered, however it does not appear that sub- 
stitution of a refractory material for Plexiglas, without 
a change in mass injection geometry will completely al- 
leviate the ablation problem. The thermal diffusivity and 
ablation temperature of a refractory such as quartz or 
boron nitride are much higher than of Plexiglas, which means 
that for a given heat flux it will take the surface of the 
refractory longer to reach its ablation temperature. Con- 
sider for example a semi— infinite slab of insulator at 
initial temperature T^. At t = o let the surface be subjected 
to a flux of energy P Q W/cm 2 which thereafter remains con- 
stant in time. The surface temperature T, as a function of 
time t, is given by^ 82 ^ 


T-Tl 


2Po 

(jt k f>Cp)V»- 


where k, , and are the thermal conductivity, density, 
and specific heat of the slab. If the ablation tempera- 
ture T* is reached at t = t* one can write 


£*= t*-tJ 

Thus the ratio of the times to ablation of two different 
materials is given by: 


M = k, P , c„ (Tf- Tj f 
tl k t p*C„ t \Tz*-Ti) 
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where subscripts 1 and 2 identify the two different materials. 
Denote quartz by subscript 1 and Plexiglas by subscript 2. 
Then: 


k 1 = 7.96xl0 -3 joules/sec cm°C k 2 = 1.88xl0 -3 joules/sec cm°C 
3 

jOj = 2.66 g/cm 


Pa 


= 1.2 g/cm' 


Cp- L = 1.01 joules/g C 
T x * = 1600 °C (78) 


Cp 2 = 1*47 joules/g C 
T 2 * = 460 °c (83) 


Thus 


+* 


« 80 


for T^ = 25 °C. Kerr-cell photographs similar to those in Fig.A-1 
indicate that ablation of the Plexiglas near the base of the 
cathode begins about 5-10 psec into the current pulse. Thus, 
if quartz were used instead of Plexiglas, it would begin to 
ablate there sometime between 400 and 800 psec after discharge 
initiation, all other conditions being equal. For the remain- 
der of the current pulse the ablating surface area of the quartz 
would continue to grow, albeit more slowly than for the Plexi- 
glas. It is not known, however, how the current density ad- 
jacent to the quartz insulator, and thus the rate of energy 
deposition to its surface would be affected by the lack of 
either argon or ablated material between the injectors during 
the transient heating time. 


A-4 TRANSIENT EFFECTS ON SPECIES STRUCTURE 

Since the species structure of the exhaust is a result 
of the ablation of the Plexiglas insulator, it follows that 
during the first 150 psec of operation the structure will re- 
flect the growth of the ablation rate. The Kerr-cell photo- 
graphs of Fig. A-l confirm this. At 40 psec, just after sta- 
bilization of the current pattern, no structure is yet visible 
in the plume downstream of the anode or in the conical luminous 
region between cathode and anode. By 100 psec, however, the 
familiar jet structure is clearly visible, and at 150 psec it 
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exhibits its final steady configuration. The increase of 
ablated flow rate during this time effectively lowers the 
average molecular weight of the exhaust (the average mole- 
cular weight of the ablation products is about 6.6 as com- 
pared to 40 for argon) , and, it is believed, leads to an 
increase in the average exhaust velocity. The latter can be 
justified on the basis of experimental results of Malliaris ^ 
who found that regardless of the propellant used, the 
kinetic energy of the exhaust of an MPD accelerator re- 
mains constant, so that the average exhaust velocity scales 
inversely as the square root of the average molecular weight. 
The increase in the average exhaust velocity should be de- 
tectable as an increase in the back emf of the arc. Terminal 
voltage measurements. Fig. A-5, confirm this. Figure A-5c, 
obtained at 16 kA and 6 g/sec shows a monotonic rise of about 
32V during the first 150 nsec, after the drop from bank vol- 
tage (4KV) . The other voltage signatures, obtained at dif- 
ferent operating conditions, indicate that more pronounced 
voltage increases occur under starved conditions while flat- 
ter signatures result under overfed conditions. The observed 
transient voltage change at any condition is, however, not 
entirely due to the increase in back emf. Part is believed 
to correspond to the energy the discharge must expend in 
pyrolyzing and dissociating the material it ablates. For 
example, indications are that for the 16 kA, 6 g/sec case 
the rate of ablation is approximately 1 g/sec. ^ ^ For Plexi- 
glas the energy of surface depolymerization and gasification 
is about 27 kcal/mole, while the energy of decomposition 

is 1470 kcal. ^ Thus an ablation rate of 1 g/sec requires 
a power expenditure of about 6 x 10 4 W, which at a current of 
16 kA, corresponds to a voltage of approximately 4V, which 
is about 12% of the increase and only about 2% of the total 
steady state discharge voltage. Most of the voltage change 
during the first 150 |isec would thus appear to be a result of 
increased average exhaust velocity. Clearly, direct measure- 
ments of the ablation rate are necessary for a more accurate 
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assessement of the relative contributions of back emf and 
insulator ablation to the terminal voltage. 
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APPENDIX B 

THEORY OF THE FABRY-PEROT INTERFEROMETER 
B— 1 INTRODUCTION 

The Fabry-Perot interferometer, shown schematically in 
Fig. B-l , consists essentially of two transparent glass or 
quartz plates^ whose faces are in principle perfectly plane 
and parallel. Each of these faces is covered by a reflecting 
coating of reflection coefficient R, transmission coefficient 
T, and coefficient of absorption A = 1 - R - T. If the plates 
are illuminated by a convergent (or divergent) beam of mono- 
chromatic light the multiple reflections off the inner plate 
surfaces give rise to a system of interference fringes at in- 
finity which when focused on a screen by a lens appear as a 
set of concentric luminous rings. 

A ray of light incident on the plates at an angle 0 , 
gives rise through successive reflections to an infinite num- 
ber of parallel transmitted rays of decreasing amplitude. The 

path difference between any two successive transmitted rays 
is (24,41) 

S = 2 pd cos© 
and thus the phase difference is: 

<t> = 2 rr i- = 2 7 Tp 
A 

where p is the refractive index of the gap between the 
plates, d is the plate separation, X is the wavelength of the 
incident light, and p= ( 2pd cos&)/\ is the order of inter- 
ference. If p is an integer, i .e . ; the path difference 
between successive rays is an integral number of wavelengths, 
the transmitted rays interfere constructively and the incident 
ray is transmitted. If this condition is not exactly ful- 
filled the multiple rays interfere destructively, the con- 
structive interference of a great number of vibrations being 
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very critical, and the incident energy is reflected back to 
the source. 

If E q is the amplitude of the incident light the complex 
amplitudes of the transmitted rays are 

E q T, E q TR e -1 ^ , ...E Q TR n e -in ^ . . . 

which form a geometric progression with common factor R . 

The sum of these amplitudes at the focal plane of the imaging 
lens is therefore 

£ = - E * r , 

T t-Re* 

The intensity Irp = [E^j °f the transmitted light is thus 

LT Z 

|l - 

2 

where I q = |E q | is the intensity of the incident beam. 

The intensity thus varies with <f> according to the 
periodic function: 

A ' ' < > (ft? “-‘i 

(24 41 42) 

which is known as the Airy function. ' ' Hits general behaviour 

is shown in Fig. B-2, for various values of the reflectance 
R. If R is small the interference fringes are broad and 
indistinct, whereas if R is close to 1, the fringes are 
very sharp. This and other factors affecting fringe sharp- 
ness will be discussed in more detail in a later section. 

Since ^ is a function of incident angle 0 , optical 
thickness pd and wavelength A so is the function . For 
a given pd , A 1 is a function only of & and A so that a 


given by 




T-t = 
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monochromatic source of light gives rise to a unique system 
of concentric interference fringes. When the light source is 
not monochromatic the interferometer will produce a distinct 
ring system for each different wavelength present. For a 
given order of interference the diameter of the resulting 
circular fringe will be a function of the wavelength. In ef- 
fect then, the Fabry-Perot causes a spatial separation or dis- 
persion of wavelengths. 

B-2 DETECTION METHODS 

Clearly, the simplest method of recording the interfer- 
ence fringes is to place a photographic plate at the focal 
plane of the fringe imaging lens. This method is generally 
used with interferometers of fixed plate separation, better 
known as etalons . A more versatile detection system which 
lends itself readily to time-resolved spectroscopy, is one 
which uses photoelectric detection, with the interferometer 
operating in the so-called scanning mode . In 

this technique the photographic plate is replaced by a dia- 
phragm which isolates a region of points which receive rays 
of the same wavelength. In principle this diaphragm should 
be an annular slit, however a pinhole located at the center 
of the fringe system is generally employed. The pinhole can 
be considered as an annulus of zero inner radius. The inter- 
ferometer and diaphragm together act as a filter or mono- 
chromator. To make a spectrometer a photoelectric detector, 
coupled to a recorder, is placed behind the isolating dia- 
phragm, as shown in Fig. B-3. The transmitted wavelength 
at the pinhole is varied continuously by varying the optical 
thickness p.d of the interferometer. One way of doing this 

is to vary the plate separation d by means of a piezoelec- 

(44 45) 

trie transducer applied to one of the plates. ' 

If d is made to decrease, the fringe system contracts, each 
fringe decreasing in radius until it disappears at the center. 
If several discrete wavelengths are present their separate 
ring systems are thus swept by the pinhole successively, so 
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that in effect the spectrum is scanned. 

At the pinhole, which selects collimated light normal 
to the plates, the wavelength of maximum transmission is a 
linear function of the plate separation. However, for any 
given wavelength the device has an infinite number of trans- 
mission modes such that for each amount A/2 that the plate 
separation changes, the given wavelength will again be trans- 
mitted (Fig. B-4) so that the transmitted intensity function 
as recorded by the detector will consist of an infinite num- 
ber of equally spaced fringes, i . e : the Airy function. 

Conversely, if the wavelength A at which transmission 
occurs is changed at fixed plate spacing, the instrument 
will transmit again when A has changed an amount 
^Af = A 2 /2|id. (24 ' 41) This is equivalent to the ring system 
changing size as A changes until the p'th ring at wave- 
length A + aA f overlaps the (p+l)th ring of wavelength A . 
This wavelength change is called the free spectral range. 

Any two wavelengths separated by an integral multiple of 
4A f will be transmitted at the same plate separation. 

(Fig. B-4) Consequently the wavelength range that may be 
examined unambiguously with the Fabry-Perot interferometer 
is limited to aA f. (In terms of wavenumber or = 1/A , 
the free spectral range is given simply by A<T f = l/2pd) . 

B-3 FACTORS INFLUENCING THE FABRY-PEROT SPECTROMETER 

By definition, the instrumental function of a spectro- 
meter is the curve which would be recorded if the instrument 
were illuminated by a perfectly monochromatic source. An 
infinitely thin spectral line is always rendered as one with 
finite width by any real spectrometer. The overall instru- 
mental broadening function is a convolution of the various 
individual broadening functions present in the device. ^ 42 ^ 

In a scanning Fabry-Perot spectrometer there are three 
distinct contributions to line broadening. In the first place, 
even if the instrument is otherwise perfect the interference 
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fringes will have a limiting finesse, that is, a finite width. 
The finesse is defined as the ratio of fringe separation to 
the width at half intensity of a fringe. This theoretical 
finesse depends only on the reflectivity R of the plates 
and is called the reflectivity finesse, F . In reality the 
effective finesse of the fringes is also limited by the un- 
avoidable imperfections on the plate surfaces which cause the 
plate separation to not be strictly constant across the di- 
ameter of the plates. Furthermore, the necessarily finite 
size of the pinhole also causes a decrease of the finesse. 

The three effects just described can be characterized by three 
partial instrumental functions, A^, A^ and in that order. 
The resulting total instrumental function is a convolution of 
these three partial functions: A T = a 1 * A 2 * A 3 . 

Effect of Plate Reflectivity 

A Fabry-Perot interferometer with perfectly flat plates, 
illuminated by monochromatic radiation of wavenumber a has 
a transmission function 



where ss 


- ___ 

(l~R) z 

T a 

PW 


as was shown previously. If the pinhole scanning aperture 
is infinitely small, the scan of the interferometer will 
register the above function. The ideal instrumental function 
A l is thus the Airy function, as shown in Fig. B-5, where it is 
plotted against wavenumber, 0 - . it consists of a series of 
maxima whose width at half height is and which are 

separated by the interval The latter is the free spectral 

range which in terms of the wavenumber is l/2pd. . The half- 
width of the Airy fringes for values of R not less than about 
0.8 is well approximated by (24,41, 42) 


&l •— 


l-R 

n\/7C 
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The ratio ~ F R is ^e reflectivity limited 

finesse of an otherwise ideal instrument. It increases 
rapidly as R approaches 1. With plates of very high re- 
flectivity a finesse of several hundred can in principle be 
attained. The maximum transmission factor is theoretically 
independent of R if the absorption coefficient A of the 
reflecting coatings is negligible: 


i . e : 




/ 


if A=1-T-R = 0 


However, if A / 0, will decrease rapidly as R approaches 

1. (Fig. B-6) It is desirable thus to use reflective coatings 
with as low an absorption coefficient as possible. Typical 
values obtained in practice with multilayer dielectric films 
lie in the range 0.05 - 0.2%. 


It would seem, on the basis of the preceding arguments, 
that as high a reflectivity as possible should be used to maxi- 
mize the instrumental finesse and thus its resolving power. 
However, the true finesse is always ■ limited by departures from 
perfect surface flatness of the plates, and the reflectivity 
finesse F R should not be greater than a certain value, which 
is determined by the quality of the plates. 

Influence of Surface Flatness 


To evaluate the reflectivity limited instrumental function 
A 1 (<r ) the interferometer plates were assumed to be perfectly 
plane and parallel, so that the plate separation was everywhere 
exactly equal to d. In actual fact however, it is impossible 
to manufacture plate surfaces which are flat enough that sur- 
face defects may be neglected. This is easily verified ex- 
perimentally by illuminating the interferometer with a normal 
collimated monochromatic beam of light and adjusting the plate 
separation for maximum transmission. On a screen placed im- 
mediately behind the interferometer a uniform circle of light, 
with a diameter equal to that of the incoming beam or of the 
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clear plate aperture, whichever is smaller, should be seen 
if the plate separation is constant everywhere. However, in 
reality the circle of light will not be uniformly illuminated 
but will consist, in general, of irregular regions of light 
and darkness such as shown in Fig. B-7, which is observed with 
the Fabry-Perot interferometer used in the experiments de- 
scribed in Chapter 4. This proves that only over certain 
areas of the plates is the separation exactly that required 
to transmit the monochromatic beam of given wavelength. If 
the wavelength is varied, or if the average plate separation 
is changed, the luminous zones will change shape. These 
variations characterize the surface flatness defects of the 
interferometer plates. 

The problem is usually analyzed by considering an imper- 
fect interferometer to consist of a juxtaposition of small 
interferometers of different plate separation, where the total 
transmitted light flux is the sum of the individual fluxes 
transmitted by the elementary interferometers. The exact 
instrumental function A 2 (<J) for various plate defects can 
be obtained analytically. For example, if the 
plates have a spherical surface the function is rectangular, 
if the plates suffer random micro-defects the function is 
Gaussian, and if one plate is tilted relative to the other 
the function is parabolic. In actual fact the plate sur- 
faces, and thus the plate separation, suffer from a combination 
of such defects. It can be shown, however, that regardless of 
the actual forms of the various defects, the instrumental func- 
tion has a width that is related to the mean deviation from 
(43) 

flatness. The flatness limited finesse is a simple func- 

tion of this deviation and is given by: 

P m / _ A 

Zcr&E ” ZA& 

where CkB. represents the mean deviation from flatness. For 
example, if the plates are flat to A /100, the flatness 
limited finesse F^ = 


50 . 



142 



SURFACE FLATNESS OF INTERFEROMETER PLATES 


FIGURE B-7 

AP25-4857 
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Effect of Finite Pinhole Size 


The apparatus function A^(G) characterizing the in- 
fluence of a pinhole of finite size is obviously of rectang- 
ular form. This can be understood by considering the pin- 
hole in its more general form of an annulus, and assuming 
the fringes projected on the focal plane to be negligibly 
thin compared to the annular aperture. During a scan of the 
interferometer the flux transmitted by the aperture remains 
at a constant value as long as an interference fringe is 
moving within the annular slit, but goes to zero abruptly as 
soon as the fringe moves out of the annulus. The width ^<*3 

of the rectangular function A-.(<J) can be easily related to 

J (42) 

the angular size of the pinhole: 


A 


OC o 


where o( is the angle subtended by the radius of the pinhole at 

p 

the imaging lens, i.e : o< = , where P is the pinhole di- 


2 f 

ameter and f is the focal length of the lens, 
finesse is consequently 

tr A&F I 

AG 3 

Overall Instrument Function 


The pinhole 


fJidoC 


The overall instrument function is the convolution of 
the three individual functions described above: (42) 


A = A.. * A * Ao 
T 1 2 J 

This function is similar to A^ , to the extent that it con- 
sists of a series of equidistant maxima with half-width AG^. 
and which are separated by the free spectral range AGf (Fig. 
B- 8 ). The two important characteristics of A T (<5j are the ordi- 
nate "V of the maxima and their width at half height, Ag T . The 
luminosity of the spectrometer is proportional to while 

its effective resolution 1R varies inversely as AG^ , i.e. : 

■x. ——.This resolution can in principle be increased without 

AG r 

limit. That is, the width of each of the functions A^ {G) 

A 2 (<jj A ^(G) can be made as small as desired. This is 



FIGURE B-8 
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OVERALL INSTRUMENTAL FUNCTION , A T (cr) 
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evident for A 3 which approaches zero width as the pinhole 
diameter approaches zero. It is also true for A^ and , 
because = A <T f /F r and = A S /F . Since the 

free spectral range A = l/ 2 pd . it is sufficient to in- 
crease the plate separation d to decrease A and hence 
Atf*^ and A<T 2 and consequently Ac^. In principle then, the 
effective resolution <T /Affip can be easily increased 

without limit. This is one of the primary advantages of the 
Fabry-Perot interferometer. 

A knowledge of the exact total instrumental function 
A t (cr) requires knowledge of the three component functions 
Af, A 2 , A^ • It is impossible to precisely determine the 
function A 2 due to plate surface imperfections, nevertheless 
its width A can be determined exactly. Fortunately, the 
values of Z^.and A 0 “ T depend much le^| 3 y n the exact form of 
A 2 (CT) than on its half-width A . It is known that the 
breadth of the resulting function A T (ff) is always greater 
than the largest of the widths A^, A<T 2 , ^*3 of the ex- 
ponent functions, but less than their sum. To a first ap- 
proximation one can write: 

A<5y 2 * Aff , 2 +• A cr* ■+■ A 

What remains to be done now is to choose the characteris- 
tics of the interferometer in such a way as to maximize both 
the resolution and luminosity. 

Importance of Plate Flatness 

Since the breadth A<T T of the total instrument function 
is always greater than the largest of , A and A ^ 

it remains to determine which of the latter is the most im- 
portant. It is evident that A^ can be chosen as small as 
desired by decreasing the size of the pinhole sufficiently. 
Furthermore, modern reflecting coatings made of multi-layer 
dielectric films can be made to attain reflectivities very 
close to unity without appreciable absorptivity (e . g . R « 99 . 8 %, 
A «s0.05% in the red end of the spectrum). The reflectivity 
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finesse F R thus can be made as large as desired, at least in 
the visible part of the spectrum, achieving values of the 
order of several hundred. On the other hand, it is virtually 
impossible to manufacture plates with defects smaller than 
A/200. Furthermore, the subsequent deposition of dielec- 

tric coatings cannot be made without degrading this figure to 
about A/100. As a result the flatness limited finesse 

cannot exceed a value of about 50 and consequently it is the 
function A 2 due to surface defects which is the broadest. (It 
should be noted at this point that in an instrument which is 
scanned very rapidly the moving plate may bow or vibrate, due 
to the acceleration at the beginning of the scan, to an ex- 
tent which can markedly decrease its flatness figure, thus de- 

(44) 

creasing the limiting finesse F D even further. ) 

Using the concept of finesse one can define for the spec- 
trometer system an effective total finesse, defined by 
f t = 4<r f /d <r T where is the free spectral range and 

the half-width of the function A T (ff). It follows, then, that 
the total finesse F^, is always less than the finesse F D due 
to surface defects. It is obvious that the latter has the 
role of a limiting finesse which can only be approached. 

Choice of Plate Reflectivity and Pinhole Size 

It remains to choose the widths & and that is 

the reflectivity finesse, F„, and the pinhole finesse F„ 

-K. P 

in such a way as to optimize the function A T (G) with respect 
to transmission, ZT , and half-width A (S^ compatible with the 
limiting flatness finesse F Q . 

Consider the partial function = A 1 * A 2 wh; ‘- ch is the 
function of a real interferometer, accounting for plate flat- 
ness limitations, but neglecting for the moment the pinhole 
effects. Its form is analogous to those of and A T (Fig. B-5 
and B-8) . The width ^ tf 12 of its maxima can be characterized 
by a finesse F g = ^ <r^ 2 which is simply the effective 

finesse of the interference fringes themselves. The ordinate 



147 


T e of the maxima of A^ 2 represents the maximum transparency 

of the interferometer alone ( i.e . minus the pinhole) . It can 

be shown that A 0", „ (or F ) and Z” are functions ofAs, , of 

the Airy function, i . e . are functions of the reflectivity 

finesse F R . These functions are shown in Fig. B-9, plotted 

against the ratio of the reflectivity finesse to the flatness 

limited finesse, F R /F D i '{'he actual shape of these curves 

depends on the exact shape of the flatness function A 2 (ff) ; 

the curves of Fig. B-9 are average curves which do not deviate 

significantly from those obtained from specific functions A 2 . 

The curves represent the behaviour of the ratio of effective 

finesse to flatness limited finesse, F e /F D# and the ratio 

of the maximum effective transmission Z" to the maximum 

2 2 ® 

transmission Z" - ^ = T /(1-R) of a perfect interferometer. It 
can be seen that as F /F^ increases from zero the ratio 
F e /F D increases rapidly at first, attaining a value of about 
0.7 at F r /F d = 1. From there on the shape of the curve de- 
creases rapidly and F e /F D approaches 1 asymptotically for 

large values of the ratio F D /F_. 

R D 

The maximum transmission ZT of the real instrument, on 

e 

the other hand, decreases monotonically as the reflectivity 
finesse increases. Equal to unity at F rj /F r = o, i.e. for the 
case of a perfect interferometer with A<3" 2 = 0, (F^* oo) , the 
ratio r e /r M is less than unity for all non-zero values of 
F r /F d ( i.e . of A <* 2 ) . This means that the transmission fac- 
tor of a real interferometer is always less than that of one 
with perfect plates. This is so because due to surface de- 
fects only part of the plate area transmits light, as was dis- 
cussed previously. It is evident from Fig. B-9 that there is 
nothing to be gained in using a reflectivity finesse F much 
greater than F R because the effective finesse does not in- 
crease that much while the transmission factor decreases 
rapidly. Thus a ratio of F r /F d much greater than unity will 
lead to a large loss of luminosity without appreciably im- 
proving the resolution. The optimum value of F„ is then one 

(43) R 

close to the value of F^. 
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EFFECTIVE FINESSE AND TRANSMISSION FACTOR 
AS FUNCTIONS OF FLATNESS FINESSE F D 


FIGURE B-9 

AP 25 -4859 
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A similar argument on the effects of pinhole size leads 
to much the same conclusions, i.e . the pinhole finesse should 
also be of the same magnitude as F^. Thus the angular size 
of the pinhole radius should be;( 42 ) 



Summarizing then, to obtain the largest possible lumi- 
nosity for given resolution, the magnitudes of the individual 
finesses F R F R , and F p should be chosen of equal magnitude. 
Since the quality of the plate surfaces defines F , the choice 
°f F R and F p is fixed. Thus the overall instrumental finesse 
is given by F T «0.6 F R and the effective resolution is given 
by Yl * 0 . 6 # Q where &o is the resolution obtained if the 
functions A^ and A^ have zero width while the function A 2 has 
width A<7 2 - Finally, the maximum value ZT-p of the func- 
tion AfpCo - ), which determines the luminosity of the instrument, 
is given approximately by £^*0.6 where ZT M = T 2 /(l- R) 2 is 
the maximum of the Airy function A p of a perfect inter- 
ferometer. 

Luminosity of the Fabry-Perot 

The luminosity of the scanning Fabry-Perot interfer- 
ometer is defined as the ratio of the luminous flux trans- 

mitted by the pinhole to the luminance L of the licrht 
( 42 ) 

source: -r 



The luminous flux is given by: 

§ = LZ r 5JX 

where ZT^is the total instrumental transmission factor, S is 
the useful area of the Fabry-Perot plates and fi is the solid 
angle subtended by the pinhole at the focusing lens. For an 
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optimized interferometer system: 



0.6 


T * 

(I 


whence : 


. 2 2 TT 

fL ^ TT o< — yo 

#v 0 

/ = ^SA * 3.6 


In terms of the effective resolution 


1R~ 0 . 6 : 


£ * 


r 


(TW « 


With dielectric film reflective coatings 


(1-R) 


usually lies in the range between 0.5 and 0.9 in the visible 
part of the spectrum. Assuming an average value of 0.7 leads 
to < 


This result demonstrates that the luminosity of the 
Fabry-Perot spectrometer varies inversely with its resolution, 
as in any other spectrometer, and is directly proportional 
to the useful surface area of the interferometer plates. In 
practice it is the loss of luminosity which prevents increas- 
ing the resolution of the instrument without limit. For a 
given resolution, however, a Fabry-Perot spectrometer has a 
luminosity typically two orders of magnitude higher than that 
of a grating or prism spectrometer, wherein lies its fun- 
damental advantage. 
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APPENDIX C 

RECOMBINATION OF AIII TO All 
C-l INTRODUCTION 

Under the conditions of temperature and particle num- 
ber density found within the environment of the MPD exhaust 

15-3 

plume (i.es T n < 2 eV, n ~ 10 cm ) the dominant mechanism 

6 S 

for the recombination of AIII to All is that of three-body 
electron-electron-ion collisions, rather than two-body radi- 
ative electron capture. In a three-body recombination 
event two electrons 'collide in the vicinity of an ion in such 
a manner that one of the electrons loses sufficient kinetic 
energy to be captured by one of the upper states of the ion. 
Once bound in one of these upper levels, which are in thermal 
equilibrium with the free electrons, the electron "diffuses" 
downward toward the ground ionic or atomic state. Between 
the upper states the probability is large that downward 
transitions will occur via superelastic collisions with free 
electrons. However, below a certain energy level (whose 
identity depends largely on the free electron density and 
temperature) the probability for radiative decay begins to 
dominate. When the captured electron has reached the ground 
state the recombination is considered completed. 

The relevant chemical equation for the recombination of 
AIII to All is 

A*"*" + @ -he — *- + e + 

where E is the energy released by the recombination and is 
R 

simply equal the second ionization energy of argon in this 
case. Some fraction of E_. is transferred to the free elec- 

K 

trons via superelastic electron-ion collisions and the rest 
is lost as radiation. The relative magnitudes of the two 
effects depends largely on whether or not the plasma is 
optically thin to the emitted radiation . ^ It 
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is not within the scope of this investigation to discuss 
in detail the theoretical foundations of electronic recom- 
bination, and the interested reader is referred to the ex- 
tensive literature on the subject. The treatment 
of even one-dimensional recombining flows with only one 
ionic specie present is a matter of considerable complexi- 
ty, 72 ^ as it involves not only the usual fluid dynamic 
equations of mass, momentum and energy conservation but 
also species conservation equations, the recombination 
equation, an electron energy equation and considerations 
of radiative effects. Rather than attempt a full analy- 
tical treatment of the problem, a semi-empirical approach 
based on the available experimental data will be used to 
calculate the fraction of recombination of AIII to All 
within the acceleration region, and estimate whether the 
energy released is sufficient to account for the flow ac- 
celeration. 

C-2 THE RECOMBINATION EQUATIONS 

The rate equation which determines the population n. . 

- - . . (69) 11 

of second ions is 

= + <**(*, 

where b n^ n g represents the number of recombination events 
per unit volume per unit time, and «n.n e the number of 
ionization events per unit volume per unit time. n., n. . , 
and n g are the number densities of first and second ions 
and electrons respectively. of is the ionization coefficient, 
which near equilibrium can be calculated from the Saha re- 
lation. - The quantity b is the recombination coefficient 
and is given by^^'®^ 

h> — 8.7S * lO^Z. 3 -=£>57 cm 3 /sec 

where Z is the ionic charge, and T g is the electron tempera- 
ture in eV. 


4%u 

71 
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At the temperatures present in the exhaust plume 
downstream of x =5 cm ( i . e . T g ^1.7 eV) small changes 
in T e can drastically affect the value of e* , the ioniza- 
tion coefficient, because of its exponential dependence 
on temperature: ©<'■'■' e 1 /kT , where E. is the ioniza- 

/ /■ a\ ® 1 

tion energy. ' For example, a decrease in T g from 

1.7 eV to 1.5 eV for = 27.6 eV decreases cx by nearly 
an order of magnitude. Thus, as the flow proceeds down- 
stream of x =5 cm the ionization rate of All to AIII 
quickly becomes very small due to the decrease in T g , 
and the ionization term in the rate equation can be 
dropped. The recombination term will also decrease with 
increasing x because of its dependence on the particle 
densities which decrease in the downstream direction; how- 
ever, the recombination coefficient, b, is temperature sen- 
sitive ( b~/T~ 9 / 2 ) so that the decrease in recombination 
rate due to the drop in density of the expanding plasma is 
compensated to an appreciable degree by an increase in the 
recombination rate coefficient due to cooling. 

The rate equation with the ionization term dropped sim- 
plifies to: 

= - i % 

= -6 

27 3 0 q / ^ 

where B = 8.75 x 10 crn ev / 2 sec . Since each recom- 

bination event results in the disappearance of one electron 
for each doubly ionized atom recombined, the rate of change 
of electron population is equal to the rate of change of 
second ion population, i . e : djJk. = d^Ai 

at at 

= - 3 



'flu -nf 

T* 


Halt 

T 4 



Therefore, 



In addition to this rate equation, the electron 
conservation and the global mass conservation equation 
are needed to take into' account the flow expansion. Equa- 
tions of momentum and energy conservation need not be con- 
sidered because the axial temperature, velocity and dens- 
ity profiles along the argon jets are known from the ex- 
perimental measurements. 

Electron Conservation: 

Consider one— dimensional flow through the control 
volume shown below, in which there is a source (or sink) 
of electrons: 



If electrons are created (or destroyed) at a rate n g , then 
the rate of electron mass addition (or depletion) within 
the volume is m e n e AAx = f> e A Ax, where m e is the elec- 
tron mass, A the flow area, and the mass density of the 
electron gas. The rate of increase of electron mass in the 

control volume is equal to the net inflow through conserva- 

( 56 ) 

tion plus the source term. Thus, 

j- t (p> A ) + j[(e* uA ) = fi A 


For steady flow — Q and therefore 


jx(f* U/ “ P> A 

Dividing through by the electron mass, m , one obtains 
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For steady one-dimensional flow the global continuity 
equation is simply 

&(r*) - ° 

where p= n^nu; n Q being the total number density of heavy 
particles and itk their mass. 


Now define an ionization fraction: f = n /n . 

e o 


Because 

all the argon atoms are either singly or doubly ionized, f 
will in general be greater than unity. Substitution of this 
expression in the electron conservation equation yields: 


J_ J-tfU'MiAu.) = A d&e. 
ix ' ' Jit 


where the left hand side has been multiplied by 
Since J - n_nK, this equation becomes 

puA At + JL A (pa A) - A dzu 


m . 
1 

m^ 


= 1 , 


The second term on the left vanishes as a result of the glo- 

dn e 

bal continuity equation. Substituting for from the 

recombination equation, and using the fact that n g = n Q f, 
so that n^ = (f-l)n Q , one obtains, finally, 

¥ = ~ B ¥ f TV^ 

Jx u 7I 9 (* 


A change of variables, x = x^, where x, is a loca- 

tion where f is known, results in: 

At = _ a X, 71? fYf-l) 
b u 77% 

which can be integrated to yield the following expression: 

*L. ds + ^ 


f + ^ Pt 1 ) f + K 

j, 
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where K t 


1 

f, 


-h Jb> t- 



and is evaluated at ^ = 1. 


The axial profiles of n , T and u are available from 

e e 

the experiments, and the right hand side of this equation 
can be integrated numerically. However, for the purposes 
at hand it is more convenient to approximate the experi- 
mental curves by simple power laws, i. e ; 

5 * 

u. = u, 

T, - X, S c 

Upon substitution of these relations in the equation above, 
and integration of the right hand side, the following re- 
sult is obtained: 


where 


- St 

C, - ~ 8.75 x I 0 ~ zl 

66 / 


S =ft -/ 


and 


5 at 2a - b 


This equation for f is transcendental and must be solved 
graphically o 

The experimentally determined flow velocity profile in 
the argon jets, shown in Fig. 4-18 can be approximated 
reasonably well by u ~ x^ in the interval from x = 5 cm 
to x = 15 cm, where the bulk of the acceleration outside the 
chamber occurs. Thus in u = u. ^ ^ b takes on the value 

, 1 s 

1/2. It has already been argued that T q ~ x is a good 

representation of the electron temperature profile; thus in 

T Q = T q1 ^ c , c takes on the value -1/2. The profile of 

heavy particle density n Q presents more difficulties. The 

contours of n in Fig. 5-2, from which n can be estimated, 

^ o 

are not very accurate downstream of x = 5 cm, and, actually, 
no contours exist for locations in the argon jet beyond 
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about x = 10 cm. In the interval 5 < x <10 cm it ap- 
pears that x which is unusual, considering that the 

lateral dimensions of the argon jets appear to increase 
linearly with x and the velocity increases as x^. 

Based on the continuity equation one would expect 

-2k 

x , not accounting for recombination. The same 
behaviour is expected of the heavy particle number density* 

-2k 

n Q ' x , although no assumptions about recom- 

bination need be made in this case. These contradictory 
indications for the functional relationship between n and 
x are likely a result of large experimental errors in the 
electron density contours downstream of x = 5 cm. In the 

, _ol^ 

present case it will be assumed that n ~ x 2 is the cor- 

o 

rect functional dependence for 5 < x < 15 cm and conse- 
quently the value of a in n = n . % a takes on the value 

q° °1 

-2k° Hence, s = 2a - b - c = -3.25. 

The conditions at ^ = 1 (x = 5 cm) are: 

n . = 1.45 x 10 15 cm -3 

ol 5 

u^ = 8.5 x 10 cm/sec 

T , = 1.7 eV 
el 

Z =2 

„ „ . a 15 -3 

n = 2.2 x 10 cm 
e 

f x =1.54 


Thus , 


j- ■+ “ 3.1 x I0' 2 ( f* is - / ) - OAO 


The graphical solution for | = 3 (x = 15 cm) yields f=1.5 
The fraction ( f — 1 ) of second ions thus drops from 0.54 at 
x = 5 cm to 0.50 at x = 15 cm. In other words about 1% of 
the second ions recombine in this axial interval. If the 
entire energy of recombination were recovered as random 
thermal energy and subsequently as streaming energy, only 
about 1 eV per heavy particle would be recovered - an in- 
significant amount. Actually under the conditions prevalent 
in the exhaust plume, more than half the recombination energy 



158 


would be lost in radiation, as can be ascertained through 

(73) 

calculations similar to those of Fraser et al. (It is 

interesting to note here that because of the dependence of 
the recombination coefficient on the cube of the ionic 
charge, the predicted recombination rate of All to AI is 
nearly an order of magnitude less than that of AIII to All, 
which is consistent with the spectroscopic observations.) 

As a check on the sensitivity of the recombination 

calculations on the exact profiles of electron temperature 

and heavy particle density, two other sample cases were 

considered. In the first case the total ion density was 

_2t 

assumed to have the profile used previously (n Q ~ x 2 ) , 
while the electron temperature was assumed to have the 
following form: T x 1 , but with T still equal to T. at 
x = 5 cm. In the second case the T ~ x 2 profile was re- 
tained, while the heavy particle density was taken to have 
the form n Q ~ x 1 , with the same conditions at x = 5 cm. 

The calculations reveal that in the first case 14% of the 
second ions present at x = 5 cm recombine downstream, while 
in the second case 25% of the second ions recombine. The 
energy recovered in either case is* however, not very signi- 
ficant because more than half the available recombination 
energy is lost in radiation. 
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